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richez 


Les découvertes et les travaux d’une pléiade de savants 
de génie ont profondément modifié les idées sur la struc- 
ture du monde physique à l’échelle atomique et cosmique. 
Cet épanouissement des connaissances fondamentales a 
ouvert la voie à des applications dont n'osaient rêver, au 
début de ce siècle, les chercheurs les plus audacieux, 
applications qui, à leur tour, ont permis la réalisation de 


nouveaux progrès scientifiques. 


ares 


Dans les allocutions qu’ils ont prononcées à l’occasion du 
jubilé de la M.B.L.E, M. J. van der Schueren, Ministre des 
Affaires Economiques, MM. F. Leblanc et M. Hublou, respec- 
tivement Président du Conseil d'Administration et Admi- 
nistrateur-délégué de la société, ont souligné l'importance 
croissante de la recherche scientifique en mettant l'accent 
sur l’interdépendance chaque jour plus étroite de la Science 


et de la Technique. 


Fidèle à l’esprit d'avant-garde qui l’a toujours animée, la 
M.B.L.E est plus que jamais décidée à promouvoir l’acti- 
vité de ses bureaux d’études et à donner son appui à la 


recherche scientifique. 


Dans cet ordre d'idées, notre Revue publiera, comme par 
le passé, des études et des résultats expérimentaux dus à 
des chercheurs appartenant soit à des organismes scienti- 
fiques soit à la M.B.L.E. Elle s'efforcera ainsi de pour- 
suivre le but que la M.B.L.E s’ était assigné en créant sa 
Revue : aider les autres chercheurs tant académiques qu'in- 


dustriels dans leurs propres travaux. 


La M.B.L.E espère ainsi apporter sa modeste contribution 
au développement de la recherche scientifique, de plus en 


plus indispensable à Pessor de l’économie belge. 


VOORWOORD ® De Manufacture Belge de Lampes et de 
Matériel Electronique werd opgericht in 1911 en viert 
in 1961 haar vijftigjarig bestaan. Tijdens deze halve eeuw 
hebben de Natuurwetenschappen een ontzaglijke opbloei 


gekend. 


De ontdekkingen en de werken van een schaar geniale 
geleerden hebben de opvattingen gewijzigd over de struc- 
tuur van de fysische wereld op de cosmische schaal en op 


die van het atoom. 


Deze ontluiking van de fundamentele kennis heeft de weg 
geëffend naar toepassingen, waarvan de meest stoutmoedige 
vorsers van het begin van deze eeuw niet durfden dromen 
en deze toepassingen hebben op hun beurt nieuwe weten- 


schappelijke vorderingen mogelijk gemaakt. 


In de toespraken die zij hebben gehouden ter gelegenheid 
van de jubelviering van de M.B.L.E, hebben de Heren 
J. Van der Schueren, Minister van Economische Zaken, 
F. Leblanc en M. Hublou, onderscheidenlijk Voorzitter 
van de Raad van Beheer en Afgevaardigde-Beheerder van 


de vennootschap het stijgende belang belicht van het weten- 


schappelijk onderzoek en zij hebben daarbij de nadruk 


gelegd op de steeds groter wordende wisselwerking 


tussen Wetenschap en Techniek. 


De M.B.L.E is, getrouw aan de avant-gardegeest \ 7 

die haar steeds heeft bezield, meer dan ooit € 

besloten de activiteit van haar studieburelen SN © 
te bevorderen en het wetenschappelijk 


onderzoek te steunen. 


In deze geest zal ons tijdschrift, zoals 
in het verleden, studies en experimen- 
tele uitslagen publiceren van vorsers 
behorende tot wetenschappelijke orga- 

pismentor tot de M.B.L.E. Aldus zal 

verder het doel worden nagestreefd, dat de 
M.B.L.E zich bij de oprichting van haar tijdschrift 
had gesteld : de andere vorsers, zowel in de academiën, 


als in de nijverheid, helpen in hun eigen onderzoekingen. 


De M.B.L.E hoopt aldus een nederige bijdrage te 
kunnen leveren tot de uitbouw van het wetenschap- 
pelijk onderzoek, dat steeds meer onontbeerlijk 


wordt voor de welvaart van het land. 
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L'utilisation pratique du calcul matriciel 


(1° partie) 


SOMMAIRE 


Après avoir souligné le role fondamental 

que joue le calcul matriciel en Physique 
Moderne, l'auteur s'attache à montrer l'intérêt 
que présente l'utilisation de méthodes 
dolgebre matricielle pour la résolution 

de problèmes que l'ingénieur et le physicien 
rencontrent dans la pratique. 


Dans ce but, il rappelle d'abord les éléments 
de cette algèbre et dégage ensuite 

le caractère d'opérateur que l'on peut 
conférer aux matrices de correspondance. 


Enfin. il traite deux exemples très différents 
qui illustrent la mise en œuvre des méthodes 
matricielles et en montre la souplesse. 


CDU: 512.831.004.14 


INTRODUCTION 


Le calcul matriciel a été imaginé dans le but de systématiser les opérations algébriques portant sur 
des quantités groupées en complexes auxquels une certaine individualité est reconnue. 


Le type de complexe le plus simple est le vecteur, être mathématique conçu pour caractériser entiè- 
rement certaines grandeurs dirigées telles que vitesse, force, champ électrique, etc... La notation vectorielle 
apporte beaucoup de clarté et de concision à l'expression de nombreuses lois physiques, car elle permet de 
passer sous silence le choix du système de référence dont le rôle ne peut, physiquement, être tenu pour 
essentiel. 


L'introduction de la notion de tenseur, autre complexe à plusieurs composantes, est venue d'une pré- 
occupation analogue. Lorsqu'il s'agit, par exemple, d'exprimer la relation entre deux vecteurs tels que 


déplacement électrique D et champ électrique E, dans un milieu anisotrope doux, on est amené à écrire 
les trois relations suivantes : 


| m» == Exx E; =F Exy E; ae Exz E: 
p: a E Fr + Eyy Es ae Eyz ES (1) 
l D, = Ezx E; sp Ezy ES de Ezz E, 


Le contenu physique de ces relations doit évidemment être indépendant du système de référence ; 
lorsque cette exigence est formulée correctement, on est conduit a reconnaître que l'ensemble des neuf 
qUanttes ote ex rose ‚ej, Prises dans un certain ordre, se comporte, lors d'un changement d'axes, 
suivant des règles bien précises. A un tel ensemble de neuf quantités a été donné le nom de tenseur, et 
plus précisément ici, de tenseur à deux indices, parce que le rôle de chaque composante est précisé par 
deux indices. 


En fait, la notion de tenseur est beaucoup plus générale, car le nombre d'indices peut être quel- 
conque; sil vaut zéro, on a affaire aux scalaires purs ou invariants (c'est le cas de la température, par 
exemple) tandis que, lorsqu'il est égal à 1, on retrouve des grandeurs étroitement apparentées aux vecteurs. 


L'utilité de l'emploi des tenseurs va évidemment beaucoup plus loin ; par exemple, on en fait un usage 
systématique dans la théorie de la Relativité Restreinte ; là, les composantes des champs E et H, dans le 
vide, se fondent en un tenseur: le tenseur électro-magnétique, qui est à considérer dans un espace 
quadri-dimensionnel, On y rencontre aussi des quadri-vecteurs tels que : impulsion-énergie. La Relativité 
Générale fait appel à la notion de dérivée absolue d'un tenseur qui, comme son nom le suggère, est 
également une notion indépendante de tout système de référence ; sans aucunement approfondir ce point, 
il est permis d'affirmer que, dans ce dernier exemple, le concept de tenseur est quasi irremplaçable ; 
y renoncer constituerait un non-sens aux points de vue physique et mathématique. 


Enfin, signalons que la Mécanique Quantique fait appel à la notion de spineur ; celui-ci sert à carac- 
teriser les fonctions d'onde des particules ou systèmes de particules à spin. 


Comme nous l'avons dit au début, les opérations mathématiques sur tous ces complexes sont gran- 
dement facilitées par l'emploi du calcul matriciel qui, de ce fait, joue un rôle très important en Physique 
Mathématique. 


Mais à côté de ce rôle fondamental, le calcul matriciel peut également rendre des services appréciables 
lorsqu'on l'envisage en tant quoutil de travail pour la résolution de problèmes variés que l'ingénieur et 
le physicien rencontrent dans la pratique. 


On connaît l'emploi qu'en a fait G. Kron dans l'étude des circuits électriques, mettant à profit le 
caractère linéaire des relations entre courants et tensions. 


Cette dernière condition est essentielle pour que le calcul matriciel soit applicable. 


Dans la présente note, après avoir rappelé les éléments de l'algèbre des matrices, nous montrerons, 
par deux exemples très différents, que l'application du calcul matriciel à la résolution d'un problème 
donné, implique avant tout de reconnaître les groupes de quantités reliés entre eux par une opération 
linéaire ; chaque groupe de quantités est alors assimilé à un vecteur, tandis que l'opération linéaire conduit 
à une matrice. 
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I. MATRICES DE TRANSFORMATION 


Lors d'un changement de système de référence, chaque composante nouvelle, qu’elle soit d'un vecteur, 
d'un tenseur ou d'un spineur, est une combinaison linéaire des composantes initiales, Quant aux coefficients 
de ces combinaisons, ils peuvent s'inscrire dans un tableau et sont caractérisés, en dehors de leur valeur 
particulière, par deux indices qui rappellent, l’un l'identité de la composante nouvelle, l'autre celle de la 
mei initiale ; il arrive que cette identité soit elle-même spécifiée, dans le complexe, par plusieurs 
indices. 


Donc, on écrit toujours des relations du type : 


A — y D — 
u, Pel a, u, i EE ce mn (2) 


Le tableau des o: est une matrice et رہم‎ en est un élément. 
Si on opère un nouveau changement de système de référence, caractérisé par une matrice d'élément 


général ہہم‎ on écrit de façon analogue à (2), c'est-à-dire en respectant la position de l'indice variable 
dans la sommation : 
u. = e Ba io (3) 


Br a‏ سپ ey u,‏ بل 
j i‏ 


Ainsi, le passage des composantes u; aux composantes u, peut être effectué en une seule opération, 


caractérisée par une matrice d'ċlċment général yx; , défini par: 


Mey > Bi «ij (4) 


= È 


de telle sorte que : 


Uk = À 7 
j 


La matrice d'élément yx; est appelée: produit des matrices d'éléments respectifs ہم‎ et a;;. Symboli- 
quement, on écrit : 


DI A 


où les grandes lettres À, B, T désignent les tableaux correspondant respectivement aux petites lettres 
a, B, y. La formule (4) donne la règle de formation des éléments du produit de deux matrices. 


Conformément à l'usage, quand il s'agit d'écrire explicitement une matrice, on disposera les éléments 
en un tableau rectangulaire où les éléments de même premier indice figurent suivant les lignes horizontales 
et les éléments de second indice constant, suivant les colonnes; donc: 


011 ig SA Qin 
«21 AAA AOS OE Qon 
Ani ereen rennen o Ann 


Dans le produit BA de deux matrices, l'élément se trouvant à l'intersection de la ligne k et de la 
colonne j est la somme des produits des éléments de la ligne k dans la matrice B, par les éléments homo- 
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logues de la colonne j dans la matrice A, cette ligne et cette colonne étant parcourues comme indiqué 
ci-dessous : 


Une remarque importante résulte immédiatement de la règle de formation du produit: celui-ci n'est 
pas commutatif, en général. Par contre, le produit des matrices jouit de la propriété d'associativitċ. 


Les matrices du type étudié jusqu'ici sont carrées : elles comportent autant de lignes que de colonnes, 
pour la simple raison qu'il s'agit de matrices de transformation entre grandeurs comptant le même nombre 
de composantes. 


Mais, on peut généraliser : il suffit de remarquer que si, dans un produit CD, le nombre de colonnes 
de C est égal au nombre de lignes de D, on peut constituer un tableau par une règle identique à (4). 
Alors, si l'on désigne par: 


1 : le nombre de lignes de C, | 
m: le nombre de colonnes de C, 

m: le nombre de lignes de D, 

n : le nombre de colonnes de D, 


on trouve que le produit CD est une matrice à 1 lignes et n colonnes. Le schéma qui généralise (5) est 
donc le suivant : 


1 EE j 


| 
F 
>< 

| 

B 


Dès lors, il est loisible de transcrire sous forme matricielle les trois relations (1) en assimilant les ei; 
aux éléments d'une matrice carrée à trois rangées et les E,, E,, E, aux éléments d'une matrice comptant 
trois lignes et une colonne ; le résultat du produit est une matrice à trois lignes et une colonne. 

On obtient ainsi la relation suivante : 


D: Exx Exy Exz E 
D, = Eyx Eyy €yz E, (7) 
D Ezx Ezy Ezz E, 
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Les matrices : 


D, E, 
Ds et De 
D, ES 


sont appelées des matrices-colonnes ou, plus fréquemment, des vecteurs-colonnes, Dans la suite, les entités 
sur lesquelles opéreront les matrices seront des vecteurs-colonnes. 


II. VECTEURS-COLONNES 


Un vecteur-colonne u est, avant toute chose, un ensemble de n nombres u; , réels ou complexes, dont 
chacun possède une signification qui lui est particulière. On peut convenir de dire que ces n nombres 
sont les coordonnées d'un point dans un espace n-dimensionnel E, ou bien considérer ces n nombres 
comme les composantes d'un vecteur dans cet espace; de lá, vient l'habitude de parler de vecteur au 
lieu de vecteur-colonne ; cela ne présente aucun inconvénient s'il n'y a aucun risque de confusion avec 
d'autres vecteurs, par exemple des vecteurs-lignes. 


Enoncons maintenant les postulats fondamentaux qui complètent la définition des vecteurs : 
— le vecteur nul 0 est celui dont toutes les composantes u; sont nulles ; 
— deux vecteurs u) et u® sont égaux lorsque: 
gO = He و او ےک وا‎ Ai 6 im 
— la somme de deux vecteurs u® et u® est un vecteur v dont les composantes sont: 


Vi = u; 1) + u; 2) i — Il 21 00 A it 


— le produit d'un vecteur u par un scalaire a est le vecteur r dont les composantes sont : 


Ti = ab وا‎ A “en 


Il est bien clair que les deuxième et troisième postulats n'ont de sens que: 
— si les vecteurs qui y sont impliqués ont même nombre de composantes, 


— si les composantes homologues sont de même nature, ce que l'on exprime en disant que les vecteurs 
doivent appartenir au même espace vectoriel E. 


III. MATRICES DE CORRESPONDANCE 


Les matrices de transformation, avec leurs propriétés de multiplication, ont été définies plus haut 
lorsqu'ont été évoquées les règles à observer dans les opérations de changement de système de référence. 
Une première généralisation a alors été faite à partir de la règle de multiplication (4) : elle avait consisté 
à envisager des matrices rectangulaires. 


On voit mal la raison d'être de telles matrices dans les opérations de changement d'axes. Cependant, 
les matrices rectangulaires sont utiles lorsqu'il s'agit de définir une correspondance linéaire entre deux 
vecteurs appartenant à des espaces Em et E, différents. 
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En effet, toute matrice R à n lignes et m colonnes: 


Tıı T12 AO DIO ODA Tim 
Loy To ۴ ۶۷٦ Lom 
و‎ Le Ls 


permet de faire correspondre à tout vecteur u) de l'espace Em, un vecteur u™ de l'espace E, par la 
relation : 


= Rum (8) 


u) 


ou, plus précisément : 


m 
um = $ r um 702.0 ‚n (8,a) 


R définit une correspondance de Em à En. (On dit parfois aussi que R définit une application de Em 
sur E,). En général, la méme matrice R ne permet pas de définir une correspondance dans le sens inverse, 
de E, à Ex, mais une matrice de m lignes et n colonnes en définit une. 


Il peut évidemment se faire que R soit carrée [c'est le cas de la relation (7) qui établit une correspon- 


dance de Ea D] 


Dans ce cas, ona: 
v= Ru (8,b) 
ou explicitement : 


n 
O ۶۲ 1 ニル リン eee en (8,c) 
j=1 


v et u appartiennent respectivement à des espaces F, et E, , distincts ou non, mais comptant chacun 
n dimensions. 


La matrice carrée E: 


1 ad 0 
0 MAA 0 
کے رد‎ | EE ا‎ are | (9) 
0 O D ا‎ 1 
est appelée la matrice-unitċ. 
Son élément général a pour expression : 
ei; = 8; (8;; : symbole de KRONECKER) 


On a évidemment, pour n'importe quel vecteur : 
u — Eu (10) 
Enfin, pour toute matrice carrée dont le déterminant des éléments est différent de zéro, on peut définir 
une matrice inverse R telle que: 
Rot Ro = E 01 
et lon démontre que, pour la même matrice R7, on a: 


RRS EEE (12) 
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Cela permet d'écrire, à partir de (8,b), après multiplication des deux membres, par la gauche, par Rœ, 
et compte tenu de (10) et (11): 


u — Riv 


ce qui équivaut exactement à résoudre le système des n équations algébriques et linéaires figurant en (8,a). 
Remarque, — Différence entre matrices de transformation et matrices de correspondance. 


Dans un changement de système de référence, tout vecteur est considéré comme fixe et cette condition 
permet d'exprimer les nouvelles composantes de ce vecteur en fonction des anciennes. C'est ainsi qu'ap- 
paraît la matrice de transformation. 


Au contraire, lors d'une correspondance, le système de référence reste inchangé tandis que le vecteur 
subit une transformation conduisant à un nouveau vecteur, chaque composante de celui-ci s'exprimant à 
l'aide de combinaisons linéaires des composantes du premier. 


Un exemple très simple qui illustre cette différence est celui d'un vecteur, rapporté à deux axes 
rectangulaires ; 


— si lon imprime aux axes une rotation d'angle a, on a: 


Da COS o sin o U 


00 — に ピッ COS a Ws 


la matrice carrée est une matrice de transformation ; 


Fig. 1 


— si l'on imprime au vecteur u une rotation d'angle a, le résultat est le vecteur v dont les composantes 
sont données par: 


Va COS a — sin a کت‎ 


== (voir fig. 15) (14) 


We sin a cos a ti 


ici, la matrice carrée est une matrice de correspondance (ou matrice d'application). 


Dans cet exemple très particulier, il y a une parenté évidente entre les matrices carrées de (13) et de 
(14) : l'une est l'inverse de l'autre. 


En vue des applications, les notions de vecteur (vecteur-colonne) et de matrice doivent être quelque 
peu étendues. 
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IV. VECTEURS INFINIS ET MATRICES INFINIES 


Les propriétés de définition des vecteurs peuvent s'appliquer sans difficulté au cas où le nombre de 
composantes n est infini. 


Pour de tels vecteurs, les matrices se présentent sous forme de tableaux infinis du type suivant : 


Le résultat de la multiplication d'un vecteur infini u par une telle matrice peut encore être défini par 
le vecteur v de composantes : 


La somme figurant ici est maintenant une série. Si cette dernière est convergente pour toutes les 
valeurs de i, il y a un sens à parler du vecteur v, produit de R et de u. Nous supposerons cette condition 
réalisée. 


Les véritables difficultés commencent dès qu'on aborde la question du produit des matrices infinies, 
et il faut signaler que l'algèbre de ces grandeurs est particulièrement complexe (*) ; il ne peut être ques- 
tion ici d'entrer dans les détails de cette théorie. Disons simplement que l'associativité du produit et la 
possibilité de définir, pour une matrice, un inverse unique, sont liées à certaines conditions qu'il impor- 
terait de vérifier. 


V. VECTEURS-COLONNES DERIVES DE FONCTIONS D'UNE VARIABLE 
CARACTERE D'OPERATEUR DES MATRICES DE CORRESPONDANCE 


Considérons une variable indépendante z et convenons de n'envisager pour cette variable qu'une série 
de valeurs discrètes : 


CN RE N kera (15) 
n pouvant être fini ou infini. 


Par voie de conséquence, convenons également de n'envisager, pour toute fonction univoque u(z) 
de la variable z, que la série de valeurs discrètes : 


CS ee Ro موا‎ )16( 


avec 


u 16 


Dans ces conditions, si l'on assimile les u; aux composantes d'un vecteur dans un espace n-dimen- 
sionnel E, dont chaque dimension est spécifiée par l'indice correspondant i, nous pouvons dire que toute 
fonction univoque u(z) définit dans cet espace un vecteur (colonne) : 


(*) Voir par exemple l'ouvrage de COOKE and RICHARD, Infinite Matrices and Sequence Spaces (Blackie, London). 
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sa 


wg Et, wm 


Un 


Supposons maintenant que l'on effectue sur la fonction u(z) une certaine opération mathématique A 
et soit v(z) le résultat de cette opération : nous ecrirons symboliquement : 


v(z) = Aulz) (18) 
Nous n'envisagerons ici que des cpérateurs À à caractère linéaire c'est-à-dire ceux pour lesquels on a: 
Al lawl) tr ONE 2) (19) 


quelles que soient les constantes a,, a; et les fonctions u,(z), u;(z). 
En vertu de nos conventions, la relation (18) n'est à considérer que pour les valeurs de z données 

par (15), ce qui conduit au système d'équations : 
Vi = d Au ) i — 0, 1% 2, 


( aa 


Des lors, on peut se demander s'il n'existe pas une matrice de correspondance qui permette de trans- 
crire le système d'équations (20) sous forme matricielle : 


A D n (20) 


ve a ASA 20 Uo 
VA am Air WER as a u, 
e, a, ee e EE (21) 
Va ano Gti ee a Ur 
ou 
v = Au (21,a) 


A étant cette matrice de correspondance. On voit que cette dernière prend ici le caractère d'un opérateur. 
En effet, le rôle qu'elle joue par rapport au vecteur u, en le transformant en v, est analogue à celui de A, 
qui transforme la fonction u(z) en v(z). 


Lorsqu'une telle matrice À existe, nous dirons qu'elle est l'opérateur matriciel correspondant à l'opé- 
rateur À. 


Définissons maintenant la somme et le produit de deux opérateurs À et B. 
La somme $ est l'opérateur défini par: 
Sulz) = AT Buz) Au(z) E 


Quant au produit P des opérateurs A et B, il signifie qu'ayant fait d'abord agir À sur u(z) et ayant 
obtenu une certaine fonction v(z), on fait agir B sur v(z) ; le résultat est une fonction w(z) telle que : 


w(z) = B[Au(z)] 
et w(z) peut être considéré comme résultant de l'action d'un opérateur unique P sur u(z). On a donc: 
w(z) = Putz) 


Il est évident que l'ordre des opérations doit être respecté ` pour rappeler ce fait, on écrit : 


ا و ج ا 


On montre aisément que si les opérateurs A et B sont linéaires, $ et P sont aussi des opérateurs 
linéaires. 
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À partir de ces définitions, on peut démontrer les deux théorèmes suivants : 
Théorème I. — Si A et B sont deux opérateurs matriciels correspondant respectivement aux opérateurs 
A et B, l'opérateur matriciel S correspondant à l'opérateur : 


$ = ACT ا‎ 
est donné par: 


S = A + B 


Théorème II. — Si A et B sont deux opérateurs matriciels correspondant respectivement à A et B, l'opé- 
rateur matriciel P correspondant à : 


1 = ھ7 
est donné par:‏ 


pD = A 


Nous allons considérer maintenant quelques opérateurs À particuliers. 


1. Produit d’une fonction u(z) quelconque par une fonction f(z) donnée. 


a) Symbole de l'opérateur: [f]. 
b) Définition de l'opérateur : 
c) Linéarité de l'opérateur : 


[f][aui(z) + agus(z)] = afflu(z) + a,[f]u.(z) 
d) Forme des équations (20) (A = [f]): 


n (22) 


e) Opérateur matriciel F correspondant à [f]. 


Pour mettre les équations (22) sous la forme (8,c), utilisons le symbole $;; de Kronecker. En effet, 
grâce à ce symbole, on peut transcrire les équations (22) comme suit : 


Vi = > yai DA بے مت‎ al (223) 
3=0 
On voit immédiatement que la matrice F est donnée par: 
E 0 (EEE enr 0 
0 fi DRE ere ne 0 
ge ==> EE Beni ee ir A (29) 
0 0 ONE Er 7 
Lorsque la fonction f(z) se réduit à une constante c, on a: 
F = cE 


E étant la matrice-unité. 


2. Intégrale indéfinie d’une fonction u(z) quelconque. 
a) Symbole de l'opérateur: [ I 
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b) Définition de l'opérateur : 


c) Linéarité de l'opérateur : 
[ta w ala = MAINE KÉ allje) 
d) Forme des équations (20) (4 — a : 


eS a u = T ulz)edz == WO) E Sone n (24) 


e) Opérateur matriciel J, correspondant (d'une manière approchée) à 7 


Pour simplifier les calculs conduisant à la détermination de J, , nous admettrons que la suite (15) se 


` présente sous la forme : 


OP ہہ ہے‎ Ze l 2 ee, Zet eal 0 
On a donc: 
Zi =—z, ih A ee n (15,b) 
Cela étant, le théorème de la moyenne permet d'écrire les équations (24) comme suit: 
> Vo — 0 (24,a) 
2 i-1 
Vi 7707 A Eer i en (24,b) 
j=0 
avec 
ZN ET (voir fig. 2) 
En posant : 
Uj H 
uE) = u 1 (u u) 0; = a (1 F 7 8;) 
と 
avec 


Fig. 2 
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on a: 


i-1 
if > (1 Als 


[0 u; 


Dee 


の) hu; LM نی‎ EE n (24€) 


Pour mettre les équations (24,a) et (24,c) sous la forme (8,c), employons un symbole 5; défini de la 
manière suivante : 


— pour i — 0 6,0 = 0 (۰0۰27 man 
l (5; = 1 j 410,112, OTE ,i— | 
"PEA R ぷり jE ST n 


Grâce à ce symbole 8; , les équations (24,a) et (24,c) deviennent : 


EG TT ei Ope eee n (24,4) 


j=0 u; 
Ces équations montrent, qu'en toute rigueur, il n'existe pas d'opérateur matériel J, correspon- 


dant à [ | | car les facteurs entre parenthèses ne sont pas indépendants du vecteur u, 


, . . u; = Ui تہ‎ . e e と 2 D 
Toutefois, si l'on prend h suffisamment petit, le terme 一 ニーーー Q, peut être considéré comme négli- 
p p 1 P g 
Uj 
geable devant l'unité et l'on a approximativement : 
n 
Vi = 130, h (24,e) 
J=0 


Sous réserve de cette approximation, que nous admettrons toujours dans la suite, la matrice J, est 
donnée par: 


Ok” ROL ARTE 0 0 
(em میں‎ eee 0 0 
NA pe O 0 0 

jh nv (25) 
i? ۰ء و‎ NN 


3. Intégrale indéfinie double d'une fonction u(z) quelconque. 


a) Symbole de l'opérateur : [ \\ ] 


b) Définition de l'opérateur : 


e) Opérateur matriciel J, correspondant a dy ; 


li = it il 


on a immédiatement, en vertu du théorème II: 


En remarquant que : 
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c'est-à-dire : 


(dg been Es ANDER 0 0 0 
O me 0 0 0 
1 Oye ETT SER 0 0 0 
Jo = bs 2 l nee 0 0 0 (26) 
0 E 5×7٦ 1 0 0 


4. Intégrale indéfinie double du produit d'une fonction u(z) quelconque par une fonction f(z) 
donnée. 


a) Symbole de l'opérateur: | \\ f] 
b) Définition de l'opérateur : 
(ln = E dz T dz EZ) ez 
J Zo Zo 


e) Opérateur matriciel J’ correspondant a WW à 


On montre aisément que l'on a: 


c'est-à-dire : A: 
DEE E 0 0 0 
Deeg ا ا‎ EU O 0 0 0 
Es EE E are AT 0 0 0 
le. =: b Dele Ne EE 0 0 0 (27) 
ne ERE nen ix 0 0 


Apres ce rappel succinct de l'algèbre des matrices, nous allons illustrer la souplesse du calcul matri- 
ciel en résolvant deux problèmes particuliers. 


VI. APPLICATIONS DU CALCUL MATRICIEL 


V. 1 - Etude de la séparation de substances en solution par la méthode de Craig. 


Considérons deux substances À et B en solution dans un même solvant M et supposons qu'il existe 
un autre solvant N pour ces substances mais non pour M, c'est-à-dire que M et N sont non-miscibles. 


Lorsque le solvant M contenant certaines quantités de À et B est mis en contact intime avec le 
solvant N, les corps A et B quittent partiellement M et se dissolvent dans N ; quand l'équilibre physico- 
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chimique est atteint, le rapport des concentrations du corps À dans M et N est caractéristique de À, de 
sorte que : 


Cam 


G a (coefficient de partage pour le corps A) 
AN 
De même, pour le corps B: 

Csm _ ` 

Cs 


(coefficient de partage pour le corps B) 
Cette loi de partage sera supposée valable ici, sans restriction. 


Voici maintenant le schéma des opérations suivant la méthode de Craig. 


Considérons deux séries de récipients disposés suivant la figure 3. 


Au début, les deux corps A et B sont dissous dans un volume V, de solvant M, contenu dans le réci- 


pient MO; les autres récipients de la ligne M contiennent une même quantité V, du solvant pur M. Les 
récipients de la ligne N contiennent une même quantité V, du solvant pur N 


La première opération consiste à 


à effectuer les manipulations suivantes : 
a) mélanger les contenus des récipients MO et NI, 


b) laisser l'équilibre physico-chimique s'établir, 
c) séparer les liquides après décantation, 

d) les replacer dans leurs récipients respectifs, 
e) 


déplacer tous les récipients de la ligne M d'un rang vers la droite. 


Dans cette nouvelle situation, les substances À et B sont contenues dans MI et NI suivant les pro- 
portions données par : 


(1) 
Cam oi (28) 
1 ) 
AN1 
(1) 
BM1 = ۵ (29) 
G (1) 
BN1 
Wee = E Sg Des — SS یا‎ | = = 
1 EEE EL N NNee à —IH 国 
0 1 2 0 TA eee k-1 
Fig. 3 


L'opération suivante consiste à effectuer les manipulations a), b), c), d) sur les couples de récipients 
MO, NI et M1, N2 et à les faire suivre de la manipulation e). 


A l'issue de cette seconde opération, les substances À et B sont contenues dans les récipients ۰ 
NI et M2, N2 respectivement dans les proportions (en se limitant au corps A, l'indice A étant momen- 
tanément omis) : 


(2) (2) 
M1 = M2 
ب‎ a (28.a) 
(2) (2) 
NI Cno 


Après k-1 opérations similaires, la situation se présente comme indiqué à la figure 4. 
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Dans l'opération suivante k, la solution dans Mi est mélangée à la solution dans N(i + 1) (avec 
À. AU PRES k-1) ; après l'équilibre, décantation et séparation, les solutions sont replacées dans 
leurs récipients respectifs et l'on avance tous les récipients M d'un rang vers la droite. 


Avant cette opération, on avait: 


(k-21) 
= = a ۲ء - .۰ و‎ kii (28,b) 
(k-1) 
Ni 
Après celle-ci, on a: 
(9 
0— =a E 9۲ k (28,c) 
ت0‎ 
Ni 


Si Ton exprime qu'il y a conservation de la quantité de A, au cours de l'opération k, on trouve: 


WG Cr VC Nena VOY id, 11220 k (30) 
ou, en tenant compte de (28,b) et (28,c) : 
a Vi مع‎ iF Ve Gen == (a N Sr Vs) Cute i= IL 2, ہیں‎ DO DOW ET 


(les valeurs extrêmes de i étant conditionnċes par les limites de validité des relations 28,b et 28,c). 


Posons : | 
a NA 2 
ーーーーーーー- کے‎ A et — = À + = il (32) 
NM ラウ PNEU TER ae À に 
Il vient: 
Cues EG O E ee 1] (31,a) 


D'autre part, en appliquant les relations (28,c) et (30), respectivement pour i=1 et i = 0, et en 
tenant compte de ce que : 
(k-1) 


Mo = 0 
On obtient : 
ae e Gan (33) 


Il est aisé de voir que les équations (3l,a) et (33) peuvent étre transcrites en une seule relation 
matricielle, comme suit : 


(k-1) 
es oR O RM 78 946 0 کے‎ 
(k-1) 
SE A p ee (0) M2 
(k-1) 
(k) 
‘ N ea A i (0) 
a = 0 u Ms (34) 
(k-1) 
ہہ رس‎ EE 0 B 
(k) 
A 0 。 7929777 SM 0 


(k-1) , > fa z 
Cu étant égal à zéro. 


La matrice, qui apparaît dans cette relation, est une matrice carrée à k rangées. 
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Si, au lieu de s'arrêter à l'opération k, nous envisagions l'opération k' telle que k < k, la matrice 
carrée correspondante ne comporterait que k' rangées. 


Cependant, on peut écrire la relation du type (34), relative à cette opération k', sous la forme : 


ġa EP 97و‎ 0 wu 
: E ee ie A 0 
E EE Gre 
Se a? ا ہہ‎ O رش‎ 0 (34.a) 
A EE 0 
AA ka 7+509: ود‎ 0 


Si l'on désigne par Cy.” le vecteur à k composantes (p < k) : 


(p) 
M1 
(p) (p) 
M r- Cup \ k 
0 
0 


et par À la matrice carrée figurant dans (34) et (34,a), nous pouvons écrire successivement : 


C3) ۳ (k-1) 

Cu == A Cu 
(k) A (k-2) 

Cu = A E (34,b) 
(CO 7 (1) 

Cu SE? AED) Cu 


(1) D , N à 2 o 7 (1) 
Or, le vecteur Cy na qu'une seule composante différente de zéro, c'est la composante Cy, 


En appliquant les équations (28,c) et (30) respectivement pour i = 1 et pour k =1, i—0, on 
obtient : 


9 


= (0) 
M1 = ae Cmo 


Cmo étant la concentration initiale dans le récipient MO ; en la posant égale à 1, la dernière des équations 
(34,b) nous permet d'écrire : 


(k) 
M1 1 
(k) 
M2 0 
(35) 
ーー N AF 
(k) 
Mk 0 
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Explicitons la matrice Ak) La règle de multiplication (4) conduit immédiatement a: 


p? 0 BR 0 ۰ , 
2uA vi O Ca TT 0 
A? 2uA E DODO eee: 0 
0 As 2uA en ة7‎ 0 0 
کوٹ‎ je ie 7) 
lte AN RE e 0 
O EE 0ھ"‎ ep 
et l'on peut montrer, en passant aux puissances 3, 4, ...... ‚que la matrice AF-1 a pour structure: 
piet ORT Arn 0 O た 8 0 0 
(k-1) ui 2) ا ا‎ MR EE 0 
ASL = ? A i-1 ۱ (38) 
Ga! EE 0 
Ak) GIE E uit 
avec 
i (k-1 ) | 
Ġa چمچ‎ a 
(k-1-i)! i! 
Reprenons (35), on trouve: 
+ BC 
ye (ita 
= A (37) 
Cone = 


On voit que la distribution des concentrations parmi les récipients de M1 à Mk, après k opérations, 
est donnée par les termes successifs du développement de : 


ln En (38) 
Quant à la distribution des concentrations dans les récipients de la ligne N, on la tire de (28,c) : 
Cm 
Su = = 
a 
et, comme d'après (32) : 
S V, YY 
5 MEVA FISA 
on a: 
V 
Cri — v p Cui 
La distribution de Cx; est donc donnée par les termes successifs du développement de : 
Vi 
+ Air) 39 
va (u ) (39) 


À titre de vérification, montrons que la quantité totale de substance en solution dans les récipients 


M et N est égale à la quantité initiale V, (puisque la concentration CE Em effet: 
k k 
مت‎ VO MIM er Me 
il i=1 
= ۳٣ 
= V, 


Dans la distribution (38), il y a un terme maximum et le rang de ce terme est fonction de À, c'est- 
à-dire du coefficient de partage. 


S'il existe au départ deux corps A et B en solution et si les coefficients de partage a et f sont diffé- 
rents, les rangs des maxima de concentration seront différents, du moins après un nombre suffisant 
d'opérations. 


Une discussion élémentaire des résultats précédents, mais que nous ne ferons pas ici, montre que le 
rapport des concentrations dans un récipient quelconque : 


(k) 
AMi 


(k) 
BMi 


évolue de façon monotone lorsque k augmente. 
À titre d'exemple, la figure 5 représente les distributions des concentrations de deux corps À et B, 
obtenues pour k = 10, dans l'hypothèse où A, = 0,3 et An = 0,6, les concentrations initiales Camy et 


(O) 2 2 
Vimo étant supposées égales à 1. 


Fig. 5 
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VI.2 - Recherche d'une solution particulière d’une équation différentielle linéaire. 


L'équation différentielle (homogène) de la forme : 


(40) 0 7 8۹ ۶ھ 


gouverne les trajectoires paraxiales (s'écartant peu de l'axe) des électrons dans une lentille magnétique, 
c'est-à-dire dans un champ magnétique présentant une symétrie de révolution autour d'un axe rectiligne. 
z est une abscisse comptċe le long de cet axe, r est la distance à laquelle se trouve l'électron par rapport 


à cet axe; quant à f(z), c'est une fonction proportionnelle au carré du champ sur l'axe, à l'abscisse z, qui 
fait aussi intervenir l'énergie de l'ċlectron. 


Très généralement, f(z) ne diffère de zéro que dans un intervalle limité (z,, z,) ; en dehors de ce 
dernier, la trajectoire r(z) est évidemment une droite (fig. 6). On se propose de rechercher par voie 
numérique la trajectoire correspondant aux conditions initiales ci-après : 


( r(z,) = Lo 
(3262) = ës 
f (z) 
ka 
2 
Yo L, 
h h 
E Za 
Zo Zu Z, Zi る 
Th 
Fig. 6 
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En vue de tenir compte, dès le début, des conditions initiales, il y a un intérêt évident à intégrer deux 
fois l'équation (40), à partir de l'abscisse z。: ce qui conduit a: 


r (z) — ro — r。(z 一 z。) + l dz' L (2 ee ) CE = 0 (41) 


Dans les méthodes habituelles, on part d'une fonction r'”(z) dite solution d'ordre zéro et, par voie 
numérique, on calcule à l'aide de (41), où l'on remplace r(z) sous le signe dintegrale par r''(z), une 
nouvelle fonction r®(z), dite solution d'ordre un. A l'aide de cette dernière, on passe de la même façon 
à une solution d'ordre deux, puis d'ordre trois et ainsi de suite; on s'arrête lorsque la solution d'ordre p 
ne diffère plus sensiblement de la solution d'ordre p-1. 


Un tel procédé est dit itératif : il ne conduit à une solution qu'à condition d'être convergent, ce dont 
on n'est pas toujours assuré car cela dépend de la solution d'ordre zéro, r(” (z), choisie. 


Le calcul matriciel offre la possibilité d'éviter ces iterations et d'aboutir directement à la solution 
cherchée, du moins dans les limites d'erreur d'un calcul numérique qui ne peut envisager pour la variable 
indépendante z que des valeurs discrètes. 


Nous allons montrer, sur l'exemple choisi, la manière dont il y a lieu d'ordonner les opérations : 


— Adopter pour la variable z la série des valeurs : 


i lime A = سے‎ K dz' に f (2) r (z”) dz” EON on (4la) 


Zo 


— Transcrire (41,a) en utilisant les symboles A des opérateurs tels qu'étudiés en V : 


= 2 


y= Ra سے‎ o il = — din, NA او او‎ eg n (41,b) 


— Transcrire (41,b) sous forme matricielle en remplaçant les opérateurs A par les opérateurs matriciels 
À correspondants (toujours sous réserve de l'approximation en ce qui concerne J,) : 


los fo 0 OEE 0 0 0 tà 
Pen toh 0 GALA Fine A 0 0 0 Ti 
r, — ro — r。2h fo CAN و کر موا‎ a 0 0 0 i 
Ty == Me = E DN = h? DIR LA N a l-ba 0 0 0 ie 
i, == e 776 (n-1)£, (ni ee مہ‎ 0 0 In 


٦‏ ى6 
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_— Mettre (41.c) sous 


Tn 


forme explicite : 


roh 
re 2h x: خ٣‎ 
3h ہے‎ h? (Bese, eE Epi 


ro nh — h? KEIN RA E AEE HF aE 


Les calculs s'organisent donc très simplement ; la trajectoire s'obtient pas à pas, sans aucune itération 
(à l'approximation propre au calcul des differences). 


Notons enfin que les seules opérations effectives, auxquelles conduit le calcul matriciel, consistent, 
comme le montre l'exemple traité, en des additions et des multiplications ordonnées suivant des schémas 
très simples. Elles peuvent donc être facilement programmées pour être effectuées sur des calculatrices 


digitales. 
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SOMMAIRE Dans la présente note, l'auteur expose, 
dans ses grandes lignes, l'incidence 
du caractère aléatoire des phénomènes 
intervenant lors des mesures de radiations 
sur les résultats de ces mesures. 


Il examine d'abord les origines de ce 
caractère aléatoire. 


Il donne ensuite la loi de distribution 

des probabilités, relative au nombre 
d'impulsions enregistrées. Il rappelle, à ce propos. 
les éléments fondamentaux de Statistique 
mathématique, qui fournissent la justification 
des méthodes pratiques utilisées et permettent 
d'évaluer l'erreur statistique sur un 
dénombrement. 


Une brève analyse de l'influence du 
mouvement propre de l'installation de 
détection sur la précision des mesures et 
quelques considérations sur les corrections 
dues au temps de résolution de cette 
installation terminent l'étude. 
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C. — STATISTIQUES DE COMPTAGE 


I. Aspects statistiques des mesures de radia- 
tions. 


I. 1. Origines du caractère aléatoire des 
mesures. 


Une mesure d'activité d'un échantillon ou d'in- 
tensité d'un rayonnement consiste, en principe, à 
dénombrer les impulsions qu'un détecteur sensible, 
exposé à l'action des radiations étudiées, enregistre 
pendant un intervalle de temps déterminé. 


L'expérience montre qu'en répétant une telle 
mesure, dans les mêmes conditions expérimentales, 
le nombre d'impulsions enregistrées varie d'une 
observation à l'autre. 


Du point de vue de la Statistique mathématique, 
ce nombre peut donc être considéré comme une 
variable aléatoire. 


Ce caractère aléatoire est dû à plusieurs causes : 


— D'une part, le phénomène de désintégration, 
lui-même, se produit, dans le temps, sous forme 
d'événements discrets et les radiations émises 
peuvent se distribuer au hasard dans l'espace 
quant à leur direction. 


— D'autre part, les interactions de la matière 
sensible du détecteur et des radiations qui 
l'atteignent, de même que les phénomènes ulté- 
rieurs dont le détecteur est le siège, sont éga- 
lement de nature statistique. 


Nous allons examiner plus en détail ces trois 
causes. 


1°) La loi de désintégration des atomes d’une 
source radioactive. 


La loi de désintégration des atomes d’un élément 
radioactif peut être formulée comme suit : 


Pour chaque atome d'un élément radioactif, la 
probabilité dP de se désintégrer pendant un inter- 
valle de temps dt est proportionnelle à dt: 


dP A dr (1) 


Le coefficient de proportionnalité À est indépen- 
dant de l'âge et des facteurs physiques habituels 
de l'élément ; il ne dépend que de la nature de ce 
dernier et du type de réaction considéré. Ce coeffi- 
cient À, appelé constante radioactive de l'élément, 
a les dimensions de l'inverse d'un temps. 


D'après cette loi, le nombre de désintégrations, 
qui se manifestent pendant un intervalle. de temps 
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t —t donné, est nécessairement une grandeur 
aléatoire. Cherchons sa courbe de distribution. 


Supposons qu'il s'agisse d'atomes d'un élément 
radioactif dont la constante À soit très petite 
(grande durée de vie) et admettons, qu'à l'instant 
t, la source envisagée contienne un grand nombre 
N d'atomes. 


Dans ces conditions, si l'intervalle de temps 
t —t n'est pas trop important, le nombre (alea- 
toire) de désintégrations pendant t' —t est petit 
vis-à-vis de N, et l'on peut admettre que le nom- 
bre d'atomes (non désintégrés) de la source reste 
constant pendant toute la durée de l'intervalle 
t —t. 

Cela étant, quelle est la probabilité P(n') qu'un 


nombre n’ de désintégrations se produisent pen- 
dant l'intervalle t — t considéré 2 


Pour la trouver, divisons l'intervalle t —t en 


t'—t 
K ` 


La probabilité p' qu'une désintégration se pro- 


un grand nombre K de sous-intervalles 


duise pendant est donnée, en vertu de (1) 


et compte tenu du nombre total N d'atomes de 
la source, par la relation : 


t — t 
= AN 2 
p K (2) 
ou 
a (3) 
p K 
avec Sé 
E AN (4) 


La probabilité qu'un nombre i > 1 de désinté- 
— Y 
grations se produisent pendant vaut, d'après 
les lois de probabilités composées, (p')'. 


Il en résulte que la probabilité q qu'aucune 


(Re 
désintégration ne se produise pendant 
s'exprime par : 

oo | p' 
GE eT 
11 1—p' 


ou, très approximativement : 


q = l= (5,a) 


car on peut rendre p < 1, en prenant K suffisam- 
ment grand, 


On voit, dès lors, que la recherche de P(n') se 
ramène à la résolution du problème suivant : 


Etant donnée une épreuve qui ne peut conduire 
qu'à deux résultats dont les probabilités respec- 
tives p' et q' sont exprimées par (3) et (5,a), 
quelle est la probabilité P(n'), qu'au cours de 
K épreuves, le résultat correspondant à la proba- 
bilité p' soit obtenu n fois ? 


C'est un problème classique de la Théorie des 
Probabilités ; pour K> o, P(n’) est donné par 
la formule de Poisson : 


Il est bien entendu que cette formule cesse 
d'être valable lorsque l'intervalle de temps t’—t 
atteint les valeurs pour lesquelles le nombre d'ato- 
mes (non désintégrés) de la source ne peut plus 
être considéré comme restant constant pendant 
toute la durée de l'intervalle t —t. Nous revien- 
drons plus loin sur cette question. 


2°) L'indétermination de la direction de propa- 
gation des radiations ċmises par la source. 


La seconde cause de fluctuation des resultats 
de mesure est de nature géométrique. 


La condition nécessaire pour qu'une désintégra- 
tion puisse être détectée est que la radiation émise 
atteigne l'espace sensible du détecteur. 


Généralement, comme il n'existe pas de champ 
extérieur directif, toutes les directions sont égale- 
ment probables de sorte que la probabilité p que 
la radiation émise par la source atteigne l'espace 


sensible du détecteur est : 
(o) 


Seng 7 
p E (7) 


w étant l'angle solide sous lequel l'espace sensible 
du compteur est vu de la source radioactive uti- 


lisée (fig. 48). 


S: Source 
D : Détecteur 


Fig. 48 


3°) Les lois d'interaction entre les radiations et 
la matière, 


La troisième cause qui confère aux mesures leur 
caractère statistique est la nature aléatoire des 
interactions des radiations, atteignant l'espace sen- 
sible du compteur, et de la matière occupant cet 
espace. 


Nous nous bornerons ici, pour illustrer cet aspect 
des phénomènes, à considérer le cas particulier qui 
a déjà été analysé précédemment (*). 


Il s'agit de la détection de particules B— au 
moyen de photomultiplicateurs munis de plastiques 
scintillants. 


On montre que, dans ce cas, la probabilité p'', 
que l’arrivée d'une particule 8— dans le plastique 
scintillant provoque l'enregistrement d'une impul- 
sion, est donnée par: 


où Pm et Pam) ont les significations suivantes: 


Pm est la probabilité (compte tenu du rende- 
ment de luminescence propre du cristal, de 
son auto-absorption, de la transmission opti- 
que et du rendement de la photocathode) 
que l'arrivée d'une particule B— dans le 
plastique scintillant provoque l'émission de 
m photoelectrons (les particules B— sont 
supposées monoénergétiques) ; 


Palm) est la probabilité (compte tenu du seuil 
de discrimination adopté) que l'émission de 
m photoélectrons déclenche l'enregistrement 
d'une impulsion. 


I. 2. Nombre d’impulsions enregistrées - Dis- 
tribution des probabilités. 


Les probabilités p’, p et p'' définies plus haut 
sont indépendantes l'une de l'autre, 


Par conséquent, étant donnés une source radio- 
active et un procédé de détection, la probabilité p 
qu'une impulsion soit enregistrée pendant un 


intervalle de temps est exprimée par: 


Ph PBP NE 


(*) R. ScHAvES 一 Détection des rayons B— de faible 
énergie, Revue M.B.L.E, III, 2 (avril 1960), p. 108. 
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ou 


= (10) 


avec 
ANED D (€ =i) (11) 


n == 


De lá, on déduit que la probabilité P(n) qu'un 
nombre n d'impulsions soit enregistré pendant 


l'intervalle t — t est donnée par: 
— N 2 
Bin) qu, (12) 
n | 


Enfin, la probabilité P(n,, nə) que le nombre 
n dimpulsions enregistrées pendant l'intervalle 
t — t soit compris entre n, et nə, est: 


On démontre que l'on a évidemment : 
00 
P(O, ©) = 8 = 1 (14) 


Lors des mesures de radiations, ce que l'on 
cherche à déterminer c'est le taux moyen de 
comptage : 


R — 5 
1 
t—t ES 


Ce que l'on détermine effectivement, c'est une 
valeur de n ou au mieux une série de valeurs de n. 


Comment et avec quelle approximation peut-on 
en tirer la valeur de R2 


C'est la Statistique mathématique qui va nous 
fournir la réponse à cette question. 


I. 3. Eléments de mathémati- 


que ('). 


statistique 


Nous venons de voir qu'à la valeur n du nom- 
bre d'impulsions enregistrċes pendant l'intervalle 
de temps t'—t (variable aléatoire) correspond 
une probabilité P(n), donnée par la formule de 
Poisson (12). 


Plus généralement, on peut caractériser la dis- 
tribution d'une population (ensemble des résultats 
fournis par une mesure répétée un très grand nom- 


(*) C. E. WEATHERBURN — A first Course in Mathe- 
matical Statistics, Universitv Press, Cambridge (1947). 
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bre de fois, par exemple) par une suite de couples 
de nombres : 


(Pis Peer (Py a) یپ‎ 


Une telle distribution n'est pas d'un emploi 
commode. Aussi s'efforce-t-on de l'exprimer par 
quelques valeurs caractéristiques. Voici les deux 
plus importantes : 


— E(x): Valeur moyenne ou espérance mathé- 
matique de la variable aléatoire x : 


EGO 7 77 )16( 


n 


E(x) caractérise la tendance centrale de la 
distribution. 


— o : Ecart quadratique moyen de la variable x: 


o = var. ES (17) 


n 


ø caractérise la dispersion de la distribution. 
Quand la dispersion est faible, c'est-à-dire 
lorsque les valeurs de x sont groupées d'une 
manière serrée autour de E(x), l'écart qua- 
dratique moyen o est petit ; il augmente avec 
la dispersion. 


On démontre que, pour la distribution de 


Poisson, on a : 


70۲ )18( 


E Pin) a = mè = Va as‏ = اش 


P NN) 


(l'indice P rappelant qu'il s'agit de la distribution 
de Poisson). 


D'autre part, lorsque n devient grand (n > 20, 
par exemple), on peut montrer que l'on a sensi- 
blement : 


(du moins pour les valeurs de n ne s'ċcartant pas 


trop de n), de sorte que, pour n > 20, on peut 
remplacer la formule de Poisson par la relation : 


Ben 


1 2n 


(n > 20) 


275 (21) 


A TT 00913 90000000 0000 


P(n), p(n) 


0,10 p(n) Distribution normale 
P(n) Distribution de Poisson 
ñ= 20 

0,08 
0,06 
0,04 
0,02 
0 

5 10 15 20 25 30 35 40 

Fig. 49 


qui exprime un cas particulier de la distribution 
normale (distribution de Gauss) : 


1 20 
pla) = == € 
oV2r 


Dans la relation (21), on a: 


À titre d'exemple, les deux distributions (13) et 
(21), pour n > 20, sont représentées à la figure 49. 


On peut également montrer que, dans ce cas 
(n > 20), la probabilité P (n, , nə), donnée par la 
somme qui figure dans le second membre de (13), 


est sensiblement égale à p(n,, n。) défini comme 
suit : 


p(n) dn (22) 


De même, au lieu de (18) et (19), on peut 
écrire : 


l, p(n) n dn (23) 


J 


INS A p(n) (n—n)? dn (24) 


(l'indice N rappelant qu'il s'agit de la distribution 
normale). 
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Il est aisé de vérifier que, pour n > 20, la valeur 
de p(n), pour n  O, est pratiquement nulle 
(voir fig. 49) ; de sorte que: 


ce qui permet d'écrire (23) et (24) sous la forme: 


J 6 
Ex (n) = z p(n) n dn (23,a) 
= “SP SS 5 
"et De p(n) (n—n)? dn (24,a) 


Ces formules conduisent aux mémes valeurs 
caractéristiques : 


Fran) n (25) 


o = n (26) 


que les formules (18) et (19). 
Remarquons que l'on a: 
E oo 


— es 


p(n) dn = 1 


Jusqu'à présent, nous n'avons envisagé que le 
nombre total n d'impulsions enregistrées pendant 
l'intervalle de temps t —t. Or, ce qui nous inté- 
resse en définitive, c'est le taux moyen de comp- 
tage : 

= n 
R تس‎ )15( 


c'est-à-dire la valeur moyenne de la variable 
aléatoire : 


R = ーー (27) 


Si lon admet que l'erreur affectant la mesure 
de t' 一 t est négligeable, on a: 


R = E(R) = — (28) 
EL 
ee Ke 
E 
= 29 
R (=i (29: 
op désignant l'écart quadratique moyen de la dis- 


tribution de R. 


Considérons maintenant le cas d'une variable 
aléatoire x, somme ou différence de deux variables 
aléatoires x’ et x caractérisées respectivement par 


Ele J o 


On démontre que si : 


x = x + x” (30) 
on a 
E(x) = E(x’) + E(x”) (31) 
ew 
c= OER ae (32) 


En appliquant ces relations au cas ot un taux 
de comptage R est évalué par la différence de deux 
taux de comptage R, et Rə, on obtient, compte 
tenu de (28) et de (29), les formules suivantes : 


Ny n2‏ کے 


R — REC Re ee 


ta 


(33) 


Re = H > DS E 


er t) (=i) 


et 


Rappelons enfin que, dans une distribution nor- 
male, la probabilité que la valeur de la variable 
aléatoire x soit comprise entre E(x) —io et 


ى2۷ 


TABLEAU 1 


Limites 


Probabilité 


0,670 et 


G et 


20 et 


SIG; et 
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50 % (0,67 0: écart médian) 
68,3 % 
95,4 % 
09,790 


I. 2. Aspect expérimental du caractère aléa- 
toire des mesures. 


Comme nous l'avons dit plus haut, une même 
mesure, répétée plusieurs fois, ne peut fournir un 
résultat unique mais donne une série de valeurs 
numériques groupées autour d'une valeur moyenne 
qui mesure l'activité ou l'intensité cherchée. 


Nous avons vu que la probabilité P(n) d'obtenir 
un résultat n est donnée par la formule de Poisson. 
Cela veut dire que si l'on répète la mesure F fois, 
F étant très grand, on a la quasi-certitude d'obte- 
nir P(n)F fois le résultat n. 


Pratiquement, une mesure ne peut être répétée 
qu'un nombre limité de fois, de sorte que, dans 
une série donnée de F mesures, le résultat n peut 
apparaître fa fois, f, étant différent de P(n)F. 
Il est clair que la probabilité que f, s'approche de 
P(n)F, croît avec ۰ 


Le tableau Il, ci-après, montre les résultats obte- 
nus pour la mesure, répétée 100 fois, du nombre 
de particules détectées en 10 secondes à l'aide d'un 
compteur de Geiger, pour une source donnée 
d'activité constante. 


TABLEAU II 


Le tableau III ci-dessous et l'histogramme de la 
figure 50 montrent comment les cent résultats 
expérimentaux se répartissent entre les vingt 
valeurs numériques : 0 à 19. 


TABLEAU III 


2 
3 
4 
5 
6 
7 
8 
9 


Les cent résultats de mesure se groupent autour 
d'une valeur moyenne nioo égale à 8,71. 


À titre de comparaison, nous avons tracé, à la 
figure 50, le graphique donnant les valeurs de 


P(n).100 pour n = 8,71 ($). 


On voit que la répartition expérimentale est 
proche de la répartition théorique. 


Si nous admettons que celle-ci caractérise effec- 
tivement la distribution des résultats, l'écart qua- 
dratique moyen de cette distribution est alors, en 
vertu de (19): 


eo AM 209 )35( 


Cela étant admis, il serait intéressant de voir 
dans quelle mesure la répartition expérimentale 


vérifie, malgré la faible valeur de n (8,71), les 
prévisions du tableau 1 (distribution normale). 


(*) Pour le calcul de P(n), nous nous sommes servis des 
tables de: E.C. MOLINA 一 Boissons Exponential Binominal 
Limit — D. Van Nostrand, New-York (1942). 
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P(n), f, 


15 


10 


20 


Fig. 


Pour les limites : 


7 = ے 8700ھ >= mec‏ 


aon Of +9 1 


il faut faire la somme des f, correspondant à 


n = 6,7, 8, 9, 10, 11 ; on trouve 76 au lieu de 68,3. 
Pour les limites : 
n— 20 — 8,71 — 5,8 — 2,91 
EZ oO € 871 + 58:77 14,51 


il faut faire la somme des f, pour n = 3, 4, 5, 6, 7, 
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n' 8,721 
P(n).100 répartition 
théorique 
fn répartition 
expérimentale 
(u 
20 
10 15 20 
50 


8, 9,10, 11, 12, 13, 14; on obtient 94 au lieu Ge 
95,4. 


I. 3. Erreur statistique sur un dénombrement. 


En général, lors de la mesure de l'activité d'un 
échantillon ou de l'intensité d'un rayonnement, on 
ne fait pas une série de mesures dont on analyse 
la dispersion, mais on se contente d'une mesure 
unique qui a évidemment peu de chances de s'iden- 
tifier avec la moyenne n qui aurait été obtenue si 
une série suffisamment longue de mesures avait été 
effectuée. 


Dans ce cas, on caractérise quand même la dis- 
persion qu'une telle série présenterait, en calculant 
l'écart quadratique moyen : 


(on ==, 


van (19,a) 
à l'aide du seul résultat de mesure n disponible, et 
l'on convient de considérer l'écart ainsi obtenu 


comme étant l'erreur quadratique moyenne sur la 
mesure n, 


Cette manière de faire ne conduit évidemment 
pas à une évaluation rigoureuse de la dispersion de 
la série, mais ce résultat est quand même générale- 
ment acceptable, comme nous allons le montrer 
ci-dessous. 


En effet, les dénombrements expérimentaux 
portent sur des nombres d'impulsions enregistrées 
assez importants, de l'ordre de 1.000 impulsions, 
par exemple, et non sur de brèves séquences d'une 
dizaine d'impulsions, comme dans l'exemple étudié 
plus haut, qui a été choisi pour bien mettre en évi- 
dence la dispersion d'une série statistique. 


P(n) 


0 200 400 


Supposons qu'une très longue série de dénom- 
brements faits dans des conditions identiques four- 
nisse la valeur moyenne : 


me 000 
L'écart quadratique moyen de la série est : 
/1.000 = 31,6 
68,3 % de ses termes sont compris entre les 
limites : 
1.000°— 31,6. = 968,4 
12.0002 31,6 =" 1403156 
95,4 %, entre les limites : 
1.000 一 63,2 = 936,8 
1.000-+ 63,2 = 1.063,2 
et 99,7 %, entre les limites : 
1.000 — 94,8 = 905,2 
1.000 + 948 = 1.0948 


La courbe de la figure 51 représente l'allure de 
l'histogramme qui illustre cette série. L'ensemble 


600 800 1000 1200 n 


Fig. 51 
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des termes de la série est groupé d'une manière 
très serrée autour de la valeur moyenne n — 1.000. 


Supposons maintenant que la série de mesures 
ne soit pas effectuée et que l'on ne dispose que 
d'une seule mesure particulière n, donc d'un seul 
terme de la série. 


Pour fixer les idées, envisageons que ce terme n 
soit, par malchance, un terme de probabilité très 
faible, par exemple celui présentant un écart par 
rapport à la moyenne égal à trois fois l'écart qua- 
dratique moyen (supposé connu) de la série : 


w 1000  ) = 905,2 OU 1.094,8 


Dans ce cas, on trouverait : 


o = V 9052. = 30 
ou 

o = VIOO = 33 
au lieu de: 

o =~ 0ے‎ = + + 06 


qui est l'écart quadratique moyen vrai. 


Il y a 997 chances sur 1.000 pour que l'écart 
quadratique moyen vrai soit compris entre les deux 
valeurs 30 et 33 ainsi trouvées. 


Comme, du point de vue calcul des erreurs, seul 
l'ordre de grandeur importe, on voit que l'on arrive 
quand même à une évaluation de la dispersion 
parfaitement acceptable, réserve faite, bien enten- 
du, des cas où, par malchance extraordinaire, le 
résultat obtenu s'écarte tres fort de la moyenne 
vraie (inconnue) et conduise à une fausse éva- 
luation de l'erreur. 


Supposons encore, par exemple, que, dans le cas 
précédent, le résultat numérique trouvé soit 779, 
s'ċcartant de 5 fois l'écart quadratique moyen de 
la moyenne vraie. En partant de ce résultat, on 
trouverait un écart quadratique moyen égal à 


V779 = 28,1 


au lieu de 31,6; cette valeur est encore acceptable 
bien qu'un tel résultat n'ait qu'une probabilité 
négligeable de se présenter, comme le montre la 


courbe de la figure 50 (99,999 999 % des résultats 
étant compris entre les limites 779 et 1.221). 


En conclusion, on convient donc de dire que 
l'erreur quadratique moyenne (standard deviation), 
sur un dénombrement statistique n, est égale à la 
racine carrée du résultat numérique obtenu : 


Tout résultat numérique est transcrit accompa- 
gné de son erreur quadratique moyenne qui fixe 
sa précision statistique c'est-à-dire la dispersion 
d'une série à laquelle il aurait pu appartenir. 
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Pour des dénombrements de 100, 1.000 et 10.000 
impulsions, les erreurs quadratiques moyennes sont 
respectivement 10, 32 et 100 ; exprimées en pour- 
centage des dénombrements, elles sont de 10, 3,2 
et ا‎ 


Quand le dénombrement porte sur un nombre 
croissant d'impulsions, l'erreur quadratique moyen- 
ne absolue o croît comme la racine carrée du 
résultat tandis que l'erreur quadratique moyenne 
relative o, décroit comme l'inverse de cette racine : 


O = Vn (19,a) 
が 
jie TE 
n n vn 


On utilise parfois encore la notion d'erreur 
médiane ou erreur probable (probable deviation) 
pour caractériser la précision d'une mesure : 


5 = Go = UG ya (37) 
l'erreur médiane relative étant 
0,67 
Vn 
Remarque. Nous venons de voir que, si le nombre 
n d'événements aléatoires, observés pendant un 


intervalle de temps t —t est grand, on peut 
poser : 


D = (38) 


n = n (39) 


Dans le cas où n est le nombre d'atomes d'une 
source radioactive, désintégrés pendant l'intervalle 
de temps At = t' —t, on peut écrire (en supposant 
que n, tout en étant grand, est négligeable par 
rapport au nombre N d’atomes de la source) : 


a = i 2 (40) 
[voir relation (4)]. 


Or, n est évidemment égal à la diminution, pen- 
dant l'intervalle At, du nombre d'atomes non désin- 
tégrés ; on a donc: 


n = — AN 
d'où, en vertu de (40) : 
AN = —ANat 


ce qui conduit, après intégration, à la loi bien 
connue : 
(tt 
N = Ne d (0 
N, étant le nombre d'atomes de la source à l'instant 
initial t,. 


Cette loi exprime l'évolution de la source lors- 
qu'il s'agit d'intervalles de temps pendant lesquels 


on ne peut plus considérer le nombre d'atomes 
non désintégrés comme restant constant. 


Par contre, étant une loi de moyenne, elle n'est 
pas applicable pour des petits intervalles de temps 
(n faible) car elle masque le caractère aléatoire 
du phénomène de désintégration. 


II. Influence du mouvement propre de lins- 
tallation de détection sur la précision des 
mesures. 


Tout détecteur de radiations délivre des impul- 
sions électriques parasites, en l'absence d'échantil- 
lon ou de rayonnement à mesurer ; c'est son mou- 
vement propre (ou background). 


Ces impulsions sont dues à l’action des rayonne- 
ments parasites qui, inévitablement, atteignent 
l'espace sensible de l'instrument ; elles peuvent être 
dues au rayonnement cosmique, à la contamination 
radioactive des locaux où s'effectuent les mesures 
ou à celle des matériaux qui constituent le comp- 
teur lui-même. Ces impulsions sont distribuées au 
hasard dans le temps et leur existence donne lieu 
à l'enregistrement, par l'installation de numération, 
d'un taux de dénombrement déterminé, même en 
l'absence de source de radiations. C'est le taux du 
mouvement propre ou de comptage spontané, ou 
encore, dans la terminologie anglo-saxonne, le 
background de l'installation. 


Lors des mesures d'activités ou d'intensitċs de 
rayonnement, il faut évidemment tenir compte de 
cet effet parasite qui s'ajoute à l'effet que l'on 
désire mesurer. 


En général, on mesure, d'une part, le taux de 
comptage total R et, d'autre part, le taux de comp- 
tage spontané R.. 


Le taux de dénombrement dû à l'activité de 
l'échantillon est obtenu par différence : 
R = R—R, 


En supposant négligeables les erreurs sur les 
mesures de temps, l'erreur quadratique moyenne 


n LÉI | D 
sur R est Vv , n étant le nombre d'impulsions 
t 


enregistrées, pendant le temps t, pour la mesure de 
l'effet total [voir relation (29) ]. 

Vi 

L'erreur quadratigue moyenne sur R, est S 


n, étant le nombre d'impulsions enregistrées, pen- 
dant le temps t,, pour la détermination du mou- 
vement propre seul. 


L'erreur quadratique moyenne sur R. est donnée 
par la relation (34) : 


jil = y n EE Mo 
AA taż 


(34,a) 


Si les deux dénombrements sont effectués en des 
temps égaux t, l'erreur quadratique moyenne sur 
R. est donnée par: 


1 
co = = Vent 1 (34,b) 


et l'erreur quadratique moyenne relative est: 


1 


== 4 ae 
Or eve Vince (42) 
En pour-cent, elle vaut : 
100 
Op = 一 一 Vann (42,a) 
n—n, 


Si l'on s'exprime en termes d'erreurs probables, 
on écrit les relations précédentes comme suit : 


0, ٦ | سے‎ (43) 
E to? 
0,67 ees 
p = ھا ا کے‎ rb (43,a) 
0,67 
مو‎ = ت٠‎ Vm (44) 
n—n, 
67 
Dr = Ve CEE) 
n—n, 


En observant que n. = Rt, (44) peut se mettre 
sous la forme : 


0,67 me 
PRE — ۲ 3 1® 
ane | 7 
Re 


La relation (44,b) permet de déterminer pendant 
combien de temps il faut dénombrer les événements 
pour atteindre une précision p, donnée, pour un 


rapport des taux de comptage déterminé, les 


O 


temps de comptage t et t, étant les mêmes. 


= 2, (445) donne: 


Par exemple, pour 


pr WR ٤66 


Si l'on désire déterminer R。a 1 % près 
(p. = 0,01) et si le taux de comptage du back- 
ground est de 25 impulsions par minute, les temps 
de comptage pour les déterminations de R et R, 
doivent être de 8 heures 57 minutes. 
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Toutefois, les temps de comptage ainsi détermi- 
nés ne sont pas nécessairement les temps les plus 
courts pour effectuer une mesure donnée avec une 
précision donnée. Nous nous sommes, en effet, 
imposé des temps de comptage t et t, égaux ; le 
temps total de comptage peut être réduit, pour une 
précision donnée, en effectuant les dénombrements 
pendant des temps t et t, différents, judicieusement 
choisis. Si des mesures systématiques d'activités 
faibles doivent être effectuées, on se servira utile- 
ment d'abaques, pour la détermination des temps 
de comptage minima (*). 


III. Corrections dues au temps de résolution 
de l'installation de comptage. 


Tout dispositif de détection et de comptage pré- 
sente un temps mort 7. Ce temps mort est un 
intervalle de temps qui suit immédiatement le début 
de chaque impulsion et pendant lequel l'installation 
est dans l'impossibilité d'enregistrer un nouvel évé- 
nement soit que le détecteur ne soit pas lui-même 
en mesure de délivrer l'impulsion correspondante 
soit que le dispositif électronique d'enregistrement 
n'ait pas encore recouvré la sensibilité voulue. 


Pour ce qui est des compteurs de Geiger, le 
temps mort est de l'ordre de 10—* sec. Pendant ce 
temps, le compteur est incapable de délivrer une 
seconde impulsion. 


En ce qui concerne les compteurs proportionnels 
et les compteurs à photoscintillateur, il n'existe pas 
de temps mort au sens strict défini ci-dessus. 
Chaque événement donne lieu à une impulsion. 
Mais, ces impulsions peuvent se suivre de si près 
que les dispositifs électroniques d'enregistrement 
associés sont incapables de les distinguer, d'autant 
plus qu'elles peuvent partiellement se chevaucher. 
Pour les installations munies de ces détecteurs, il 
existe ainsi un temps mort apparent dont la durée 
ne dépend que du temps de montée de l'impulsion, 
de la limite inférieure de bande passante des ampli- 
ficateurs et du pouvoir de résolution du numéra- 
teur. Le temps de montée de l'impulsion due à un 
compteur proportionnel est de l'ordre de la micro- 
seconde tandis que celui des détecteurs à photo- 
scintillateur peut varier de quelques centièmes de 
microseconde à une microseconde suivant la nature 
du scintillateur et du photomultiplicateur utilisés. 
Dans les deux cas, les circuits électroniques asso- 
ciés peuvent être conçus de telle sorte que les pou- 
voirs de résolution obtenus soient respectivement 
de l'ordre de 3 à 4 fois ces valeurs. 


(*) Ces abagues ont été publiés par W. E. BROWNING 
— Nucleonics, 9, 3, p. 63. 
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Les numérateurs électroniques courants présen- 
tent des pouvoirs de résolution de l'ordre de 5 à 
10 microsecondes et peuvent être utilisés avec des 
compteurs de Geiger ; en fait, des installations pré- 
sentant un pouvoir séparateur de l'ordre de 100 
microsecondes suffiraient pour ce type de détec- 
teur. Pour être utilisés au mieux, les compteurs 
proportionnels et les photoscintillateurs doivent 
être associés à des numérateurs présentant un pou- 
voir de résolution de l'ordre de 0,5 a 1 micro- 
seconde et, dans les cas spéciaux, cette caractéris- 
tique peut être poussée à 0,05 microseconde pour 
les photoscintillateurs. 


En ce qui concerne l'influence du temps de réso- 
lution de l'installation de comptage sur les résultats 
des mesures, deux cas sont à considérer : 


— D'une part, le temps mort peut être constant, 
qu'un ou plusieurs événements se produisent 
pendant sa durée (fig. 52,a). 


— D'autre part, l'apparition d'un second événe- 
ment, pendant le temps mort qui succède à une 
première impulsion, peut allonger celui-ci d'un 
temps égal au temps qui s'est écoulé entre les 
deux événements (fig. 52,b). 


ig 52 


Dans l'un et l'autre cas, les événements qui se 
produisent pendant les temps morts, y compris 
l'événement initiateur, sont enregistrés comme un 
événement unique et les pertes de comptage qui en 
résultent faussent les résultats des mesures, Ces 
pertes sont d'autant plus importantes que le taux 
moyen de dénombrement est élevé. 


Dans les deux cas, les pertes de numération peu- 
vent être évaluées en fonction du taux de comptage 
et du temps mort, si ce dernier est connu. 


Dans le premier cas, les pertes peuvent s'ċvaluer 
comme suit. 


Soit R le taux de comptage observé au cours 
d'une mesure de durée t. Le temps total d'insensi- 
bilité au cours du dénombrement est 7R, Si R est 
le taux théorique moyen de comptage, qui devrait 


être réellement observé si le temps mort était nul, 
la perte moyenne par unité de temps est de RR. 
Le taux R réellement mesuré doit être corrigé, 


en y ajoutant ces pertes, pour obtenir le taux 
vrai Ri 


R = R+RR- (45) 
Pour des mesures suffisamment longues, on peut 
faire R’ = R et écrire la relation : 
R 
R = — 46 
I—Rr (40) 


qui permet de calculer le taux de comptage R' qui 
aurait été obtenu si le temps mort avait été négli- 
geable. 


Dans le second cas, l'analyse statistique des 
phénomènes conduit à la relation de correction : 


٢ (47)‏ .ھ00" 


Les relations (46) et (47) permettent le calcul 
des corrections dans chacun des cas. 


Mais bien souvent, les détecteurs et numérateurs 
ne présentent pas, d'une manière bien nette, l'un 


ou l'autre des modes de fonctionnement exposés 
ci-dessus. Ces deux modes de fonctionnement peu- 
vent se produire à la fois au cours d'une même 
mesure. De ce fait, à l'emploi des relations (46) et 
(48), on préfère souvent, pour l'établissement des 
corrections, l'utilisation d'une courbe d'étalonnage 
de l'instrument. 


Cet étalonnage peut se faire à l'aide d'une série 
d'échantillons radioactifs d'activité connue, grâce 
auxquels, tout en gardant toujours la même géo- 
métrie, on soumet le détecteur de l'installation à 
étalonner à un flux croissant de radiations. 


On trace alors un graphique (fig. 53) dans 
lequel on porte, en abscisses, les activités utilisées 
et, en ordonnċes, les nombres d'impulsions n rċel- 
lement enregistrées. Les erreurs étant négligeables 
pour les activités faibles, on extrapole linéairement 
la courbe obtenue aux activités élevées. Pour toute 
détermination expérimentale d'un résultat de 
dénombrement, le résultat corrigé n' peut être lu 
sur le diagramme en notant l'ordonnċe du point 
d'intersection B de la droite extrapolée avec la 
parallèle à l'axe des ordonnées passant par le 
point À correspondant au nombre d'impulsions n 
réellement observées. 


Activité 


Fig. 53 
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Quelques aspects de la technique 
de l'électrodéposition 


(1° partie) 


SOMMAIRE Dans cette première note, l'auteur passe en revue 
les grandeurs fondamentales telles que 
potentiel de diffusion, 
tension d'électrode, polarisation. 


Il rappelle ensuite les lois qui régissent 
la répartition du courant dans les électrolytes. 


Après une allusion aux méthodes 
analytiques, il s'attache plus particulièrement 
à l'examen des méthodes analogiques 
permettant de prévoir la distribution 

de ce courant. 


L'auteur expose enfin les modes d'action 

des différents facteurs physiques et chimiques 
susceptibles d'améliorer les dépôts métalliques 
au point de vue de leur uniformité. 
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I. PROBLEME DE L’ELECTROLYSE 


Le fait fondamental qui est à l'origine de la 
technique du dépôt électrolytique est bien connu. 


Si l'on plonge deux électrodes, sous tension 
continue, dans une solution d'un sel métallique, on 
observe, dans certaines circonstances, l'apparition 
d'un dépôt du métal sur la cathode. 

D'après la loi classique de FARADAY, la masse 
totale m de métal déposé s'exprime par: 


alt 
96.540 n 


p est la masse atomique du métal, en gr, 
n, la valence électrochimique du métal, 
I, l'intensité du courant, en ampères, 


t, le temps de passage du courent, en sec. 


La répartition de la couche métallique déposée 
sur la cathode dépend de deux facteurs essentiels : 


— la distribution du courant électrique dans l'élec- 
trolyte, 


— le rendement cathodique, c'est-à-dire la frac- 
tion de la densité de courant ayant effective- 
ment servi à déposer le métal (densité de 
de courant efficace). 


Bien que la connaissance de la répartition du 
courant soit essentielle pour comprendre l'incidence 
de divers facteurs sur le déroulement de l'électro- 
déposition, peu de travaux ont été consacrés, dans 
le passé, à l'étude de cette importante question. 


Ce n'est que récemment que des progrès consi- 
dérables ont été accomplis dans ce domaine grâce 
aux travaux de R. RousseLor (France) ['] 
qui a montré tout l'intérêt que pouvaient présenter 
les méthodes analogiques pour l'étude de l'élec- 
trolyse. 


Avant de décrire les procédés utilisés par R. 
ROUSSELOT, rappelons brièvement quelques pro- 
priétés des électrolytes. 


I. 1. Propriétés des électrolytes. 


I. 1.1. Potentiels de diffusion. 


Lorsque deux solutions aqueuses delectrolytes 
de nature ou de concentration différentes viennent 
en contact, il s'établit entre elles une différence de 
potentiel. 


Considérons, par exemple, deux solutions d'acide 
chlorhydrique de concentration différente. 


Avant leur mise en contact, les deux solutions 
sont électriquement neutres et, par conséquent, 
dans chacune d'elles, les concentrations en ions H+ 
et en ions Cl— ont une même valeur, celle-ci étant 
fonction du degré de concentration de la solution 
et de la température. 


Au moment où les deux solutions sont mises en 
contact, un phénomène de diffusion prend nais- 
sance et il y a passage d'ions H* et Cl de la 
solution concentrée vers la solution diluée. Comme 
la mobilité des ions H* est environ six fois plus 
grande que celle des ions CI, dès les premiers 
instants, la solution diluée acquiert une charge 
positive aux dépens de la solution concentrée. La 
différence de potentiel ainsi créée (elle n'est géné- 
ralement que de l'ordre de quelques mV) réduit 
la vitesse des ions H+ et accroît celle des ions 
CH. Lorsque les deux vitesses deviennent égales, 
un régime évolutif s'établit, au cours duquel les 
concentrations des deux solutions tendent à s'ċga- 
liser, en même temps que décroit le potentiel de 
diffusion. 


I. 1.2. Tensions d'électrode. 


Lorsqu'un métal est plongé dans la solution 
aqueuse d'un de ses sels, il s'établit, entre le métal 
et la solution, une différence de potentiel. 


NERNST a expliqué l'établissement de ces ten- 
sions d'ċlectrode en attribuant au métal la faculté 
de passer en solution sous forme d'ions. Cette 
faculté est caractérisée pour chaque métal par sa 
tension de dissolution électrolytique, grandeur 
indépendante de la nature de l'électrolyte mais 
fonction de la température. 


Lorsqu'un métal est mis en contact avec la solu- 
tion d'un de ses sels, deux tendances opposées se 
manifestent en ce qui concerne le mouvement des 
ions métalliques : 

— d'une part, sous l'effet de la tension de disso- 
lution électrolytique, ces ions tendent à passer 
du métal dans la solution, 

— d'autre part, les ions métalliques présents dans 
la solution tendent à passer dans le métal et 
cette tendance est d'autant plus forte que la 
concentration de l'électrolyte en ces ions est 
plus élevée. 


Sous l'effet de ces tendances, un transfert net 
d'ions s'opère dès les premiers instants. La dif- 
férence de potentiel ainsi créée à l'interface métal- 
solution (tension d’electrode) a pour effet d'arrêter 
ce transfert ; l'équilibre est ainsi établi. 


Partant de considérations tirées de la Thermo- 
dynamique, NERNST a montré qu'à l'équilibre et 
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sous certaines conditions, la tension d'ċlectrode est 
donnée par la relation : 


(1) 


où 
E: tension d'électrode, en volts (par conven- 


tion, le signe de E est celui de la charge 
ċlectrigue prise par l'ċlectrode), 


R: constante des gaz parfaits (8,313 joules / 
degré). 


T: température absolue, en °K. 


c : concentration ionique de la solution en 
ions-g/l, 


n : valence des ions, 


F: charge d'un équivalent gramme ionique 


(96 540 coulombs), 


C: constante spécifique du métal. 


Parmi les conditions de validité de la formule de 
NERNST, relevons les deux suivantes : 


— la solution doit être diluċe ; 


— l'électrode doit être réversible, c'est-à-dire telle 
que les modifications produites par le passage 
d'une certaine quantité d'électricité dans un 
sens disparaissent lorsqu'on fait passer la même 
quantité d'électricité dans le sens inverse. 


La mesure absolue des tensions d'ċlectrode n'est 
pas réalisable car, pour l'effectuer, par la méthode 
d'opposition par exemple, on doit nécessairement 
plonger dans la solution une deuxième électrode, 
et, dans ce cas, c'est la différence entre deux ten- 
sions d'ċlectrode que l'on mesure en réalité. Aussi, 
a-t-on convenu d'exprimer les tensions d'électrode 
en les rapportant à certaines électrodes de réfé- 
rence. Parmi celles-ci, citons l'ċlectrode à hydro- 
gène et l'ċlectrode au calomel (Hg. Cl). 


En fait, l'électrode de référence ne se trouve pas 
en contact direct avec la solution étudiée mais est 
plongée dans un électrolyte bien déterminé qui lui 
est propre. C'est entre ce dernier et la solution 
étudiée que l'on établit le contact en reliant les 
deux liquides au moyen d'un tube rempli d'une 
solution saturée de KCl (avec cette solution satu- 
rée de KCI, les potentiels de diffusion qui appa- 
raissent aux extrémités du tube sont pratiquement 
négligeables) (*). 


Si, dans la relation (1), on pose : 


c = 1 (solution normale) 


(*) On peut trouver la description de ces électrodes et de 
leur mode d'emploi dans les ouvrages traitant du pH. 
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on trouve: RT 
Be 6 تل‎ 6) 


nF 


On appelle tension standard ou tension normale 
cette valeur particulière de la tension d'ċlectrode. 


En la désignant par E, , on peut écrire (1) sous 
la forme: 


Ti 
E = Em RU In € (1,a) 
nF 


Le tableau I donne, pour des métaux usuels, les 
valeurs de E, mesurées par rapport à l'électrode 
au calomel. 


TABLEAU I 


Série électrochimique des métaux 
Tensions normales (exprimées en volts) de métaux 
purs par rapport à l'électrode au calomel. 


— 1.83 
Aluminium 一 1.55 | Argent 
Zinc — 1.05 | Platine 
Fer 0721 Or 

— 0.69 

— 0.58 

— 0.52 


Magnésium Cuivre 


Cadmium 
Cobalt 
Nickel 
Etain 
Plomb 


Hydrogène 


Les valeurs données dans ce tableau ne sont 
utilisables que lorsqu'il s'agit de métaux purs pré- 
sentant une surface de contact électrode-solution 
parfaitement saine. 


Nous verrons, en parlant de la corrosion, que, 
au point de vue pratique, on a été amené à recon- 
sidérer ces tensions. 


Revenons à la relation (1,a) qui permet de cal- 
culer E, connaissant E, et c. 


Les graphiques des figures 1,a et 1,b représen- 
tent cette relation dans les cas où on a respecti- 


vement E,>0 et E。 ご 0. 


Ces graphigues montrent que si, partant de 
c = 1, on diminue la concentration ionique, la ten- 
sion d'électrode E varie comme suit: 


— E,D 0: E commence par diminuer tout en 
restant positif, passe par zéro, devient 
négatif et croît en valeur absolue. 


— E, 00: E reste négatif et croît en valeur 
absolue. 


Variation de la tension d'ċlectrode en fonction de la concentration ionique. 


En électrodéposition, on se trouve généralement 
dans la région des courbes où E est négatif et où, 
par conséquent, sa valeur absolue augmente avec 
la diminution de la concentration ionique. 


I. 1.3. Polarisation. 


La tension d'ċlectrode dont nous venons de par- 
ler est une tension d'équilibre. L'expérience montre 
que, dans la plupart des cas, le passage d'un cou- 
rant électrique modifie la différence de potentiel 
entre le métal et la solution. On donne le nom de 
polarisation à l'écart que présente cette différence 
de potentiel par rapport à la tension d'équilibre. 


En ce qui concerne la polarisation de la cathode, 
elle résulte principalement : 


— de la diminution de la concentration en cations 
au voisinage immédiat de la cathode, 


— de certaines réactions physico-chimiques qui se 
déroulent à la cathode. 


a) Polarisation de concentration. La vitesse de 
décharge des cations étant proportionnelle à la 
densité de courant, la couche d'ċlectrolvte adja- 
cente à la cathode (catholyte) s'appauvrit d'autant 
plus vite que la densité de courant est plus élevée. 
Le renouvellement des cations dans le catholyte, 
tant par le transport dû au courant électrique que 
par la diffusion à partir de couches à concentration 
inchangée, ne suffit généralement pas à faire con- 
server au catholyte la même concentration que 
celle du reste de la solution. 


b) Polarisation chimique. La valeur des polarisa- 
tions observées ne peut s'expliquer uniquement par 
la diminution de la concentration du catholyte. 
Aussi admet-on que les diverses réactions physico- 
chimiques qui se déroulent à la cathode (formation 
et croissance des germes de réseau cristallin, for- 
mation des molécules de gaz à partir des atomes, 
etc...) contribuent à la polarisation cathodique. 


Ce qui intéresse le galvanochimiste, c'est de 
connaître, pour différents types d'électrolyte et 
pour des conditions données de température et 
d'agitation de la solution, la relation entre la pola- 
risation et la densité de courant. 


La cellule de HARING, représentée à la figure 2, 
fournit un moyen simple pour relever avec une 


Z 

u 
Z 
p 
Y 
/ 
Z 
u 
Z 
用 


MM 


Cellule de Haring 
Fig. 2 
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précision suffisante, la courbe de polarisation 


cherchée. 


La cellule de HARING est constituée par une 
cuve parallélipipédique remplie d'ċlectrolvte. Aux 


Si lon désigne respectivement par : 


extrémités de cette cuve sont placées l'anode et la 
cathode. Deux électrodes intermédiaires P, et P, 
en toile métallique, ceci afin qu'elles offrent une 
grande surface de contact avec la solution, parta- 
gent la cuve en trois compartiments égaux. 


Va: Vi; Vo; Ve les potentiels électriques de l'anode, des électrodes P, et P,, de la cathode ; 


000۲08 EPE. 


ea; en,e ررے‎ 2; & les polarisations (voir fig. 


les tensions d'électrode 


> 


d > 


2) de l'anode, des électrodes P, et P。 de la cathode; 


R la résistance électrique (électrolyte) de chacun des trois compartiments ; 


I l'intensité du courant traversant la cuve, 


on a: 
Ve ER en RI 
V, — E, — e”, RI 
Va — E, — e”,— RI 
d'où 
e, = (V,— Vi) (V,— V。) (E, — E) 
e. = (Vi —V;) (V: — V.) (E, — E.) 


Puisque les électrodes P, et P, présentent de 
grandes surfaces de contact, la densité de courant 
au voisinage immédiat de celles-ci est faible et, par 
conséquent, l'effet de la polarisation est négligea- 
ble. On peut donc omettre les termes e, , e”:, 


e's, e. Cela donne, pour les polarisations : 
ae = MVV) = VV e 


e 0 . 
e, et e. contenant uniquement les termes constants 
E, , Es, Ez, E. qui, comme nous le verrons plus 
loin, n'ont pas d'effet sur la répartition de courant. 


Les courbes de la figure 3 représentent les varia- 
tions de la polarisation cathodique (prise en valeur 
absolue) en fonction de la densité de courant, pour 
trois bains de cuivrage de nature différente : 
fluoborate de cuivre, sulfate de cuivre, cyanure de 
cuivre. Pour les électrolytes envisagés, la polarisa- 
tion cathodique est négative. 


L'examen de ces courbes montre que : 


— pour le fluoborate, la variation de la polarisa- 
tion est pratiquement nulle (grande vitesse de 
diffusion) ; 

— pour le sulfate, la polarisation croît lentement 
avec la densité de courant : 

— pour le cyanure, la courbe de polarisation monte 
d'abord très rapidement et puis tend vers une 
droite légèrement ascendante ; 
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E; ri ey = Ve 
Ez e, = V, 
= و‎ EN 


— l'élévation de la température et l'agitation de 
la solution ont pour effet de diminuer la pola- 
risation (accroissement de la vitesse de diffu- 
sion). 


—— Sans agitation 
---- 300 t/m 


0,6 


0,4 


0,2 


Courbes de polarisation de trois bains de cuivrage (cyanure, 
sulfate et fluoborate) montrant l'influence de la température 
et de l'agitation (R. Rousseror) 


Fig. 3 


Après ce bref rappel des propriétés des électro- 
lytes, nous allons évoquer brièvement les lois 
gouvernant la répartition du courant dans les élec- 
trolytes, lois qui ont des rapports étroits avec la 
théorie du potentiel. 


1.2. Lois régissant la distribution du courant 
dans les électrolytes. 


Un bain d'ċlectrolvse peut être assimilé au point 
de vue électrique à un milieu conducteur, en tout 


point P duquel règne un champ électrique E résul- 
tant de l'application de différences de potentiel 
entre les électrodes. Ce champ dérive d'un poten- 
tiel scalaire V(P) défini (à une constante près) 
par la formule classique : 


B ۰ (2) 


Pour arriver à exprimer, d'une manière simple, 
les lois qui gouvernent la distribution du courant, 
il faut faire quelques hypothèses au sujet de la 
conductivité de l'électrolyte : 


a) en tout point, la densité de courant est propor- 
tionnelle au champ électrique (loi d'Ohm) ; 


b) le milieu est isotrope au point de vue de la con- 
ductivite électrique ø (absence de stratifica- 
tion) ; 


c) le milieu est homogène dans toute l'étendue de 
l'electrolvte, (on néglige les variations de con- 
ductivité dans les zones se trouvant au voisi- 
nage immédiat des électrodes, le trajet des 
lignes de courant dans ces zones étant petit 
vis-à-vis du trajet total). 

Les hypothèses a) et b) permettent d'exprimer 


la densité de courant J sous la forme: 
J = Se grad V(P) (4) 


Cette relation peut également s'ċcrire : 


dV 
er BAe (4,a) 
VM ہیں‎ sa : 
désignant la dérivée du potentiel suivant la 
n 


normale n à la surface équipotentielle passant par 
le point P (le sens positif de n est dirigé vers les 
potentiels décroissants.) 

Si les phénomènes sont indépendants du temps 
(ou, du moins, assez lentement variables), le vec- 
teur J a un caractère conservatif et l'on peut écrire : 


dv] = 0 


Grâce à l'hypothèse c) et tenant compte de (4), 
on trouve immédiatement que : 


div grad V(P) = AV(P) = 0 


ou, plus explicitement : 


0V OV 0V 
Ox? + Oy” : 0z* 3 6) 


Ainsi, dans les hypothèses précédentes, V(P) 
est une fonction harmonique ; mais, cela ne suffit 
pas pour caractériser V(P) de façon univoque. Il 
faut encore préciser les conditions aux limites. 


Le bain conducteur est nécessairement limité soit 
par des surfaces isolantes soit par des surfaces 
métalliques. Dès lors, les conditions aux limites, 
auxquelles doit obéir la fonction V(P), sont basées 
sur les considérations suivantes : 


— En ce qui concerne d'abord les surfaces isolan- 
tes, il est certain qu'en tout point de celles-ci, 
le vecteur J leur est tangentiel; si n' désigne 
la normale en un point d'une surface isolante, 


on a: 
Jen 0 

ce qui, d'après (4), permet d'écrire : 
= EH (6) 
dn 


— En tout point P de l'électrolyte en contact avec 
une surface métallique, le potentiel V(P) dif- 
fère du potentiel électrique V. de cette surface 
par la somme de la tension d'électrode E. et 
de la polarisation e.; on a: 


VID = Ber 


E. est constant sur toute l'étendue d'une élec- 
trode. En ce qui concerne ee, la situation n'est 
pas la même à chaque instant de l'électrolyse. 
En effet, on peut admettre qu'avant l'applica- 
tion des tensions sur les électrodes, il y a équi- 
libre physico-chimique, de sorte que, tout au 
début de l'électrolyse, la polarisation étant 
nulle, on a: 


ASM 
La répartition du courant qui prévaut à ce 
moment-là est appelée répartition primaire. 


Après quelques instants, la polarisation appa- 
rait près des électrodes et modifie générale- 
ment la répartition du courant. Lorsque les 
polarisations sont stabilisées, la répartition de 
courant correspondante est appelée répartition 
secondaire. (*) 


(*) Cette notion est quelque peu théorique car les phé- 
nomènes de polarisation sont toujours plus ou moins instables. 
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Deux cas peuvent se présenter ici: 


1%) Si, pour chaque électrode, e. est indépen- 
dant de la densité de courant, le potentiel 
V(P), en tout point proche d'une élec- 
trode, prend une valeur constante, carac- 
téristique de cette électrode. 


S'il n'y a que deux électrodes en présence, 
la répartition secondaire de courant se 
confond, en valeur relative, avec la répar- 
tition primaire. 


2°) Si e. dépend de la densité de courant sur 
les électrodes, les conditions aux limites 
s'expriment par la relation : 


V(P) = Ve — E. — e(J) (7) 
J étant la densité de courant au point P 
considéré. 


Dans ce cas, la répartition secondaire est 
essentiellement différente de la répartition 
primaire. 


Ainsi donc, les conditions aux limites à obser- 
ver en tout point P situé au voisinage d'une 
paroi métallique peuvent être d'une des espèces 
suivantes : 


a) V(P) est constant et supposé connu par 
la relation : 


V(P) = V. — E. — e. 


dans laquelle n'interviennent que des cons- 
tantes. 


b) V(P) dépend du point P par la même rela- 
tion mais dans laquelle e. est fonction de la 
densité de courant. 


Nous allons maintenant discuter de quelques 
méthodes analytiques et analogiques de recherche 
de la répartition du courant. 


I. 3. Détermination de la répartition primaire 
du courant. 


Lorsque la polarisation est nulle ou indépendante 
de la densité de courant, la surface de Telectrolyte 
en contact avec une électrode est une équipoten- 
tielle. 


Au point de vue mathématique, le problème de 
la répartition primaire consiste à rechercher une 
solution de l'équation (5) satisfaisant à deux types 
de conditions aux limites simples : 
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dV 


12) — 0 pres des surfaces isolantes (con- 


dn 
dition de NEUMANN) 


2°) V donné et constant près de toute surface 
métallique (condition de DIRICHLET). 


1.3.1 — Méthodes analytiques. 


Ce problème est généralement très difficile à 
résoudre sauf dans les cas où le bain et les élec- 
trodes ont une configuration particulièrement 
simple. 


Le cas le plus simple est celui représenté à la 
figure 4. 


Tig. 4 


Il s'agit là d'un bain de forme parallélipipédique 
dont l'électrolyte est placé entre deux électrodes 
planes A, A; A, A, et C, C, C, C formant les 
faces terminales du parallélipipède. Il est clair que 
la répartition du courant est uniforme : les surfaces 
équipotentielles sont des plans parallèles aux faces 
terminales et les lignes de courant, des normales à 
ces plans. 


D'autre part, comme la répartition est la même 
dans toute tranche telle que a; a, a, ay c, c; GC 
(épaisseur unité), un diagramme plan suffit pour 
représenter cette répartition. Ce diagramme est 


indiqué à la figure 5. Les droites ee’ et ac sont 
respectivement une ligne équipotentielle (intersec- 
tion d'une surface équipotentielle du bain avec le 
plan a; a; c, c; de la figure) et une ligne de 
courant. 


Fig 


Généralement, lorsque la configuration du bain 
et des électrodes ne s'écarte pas trop d'une confi- 
guration cylindrique, le problème étudié peut se 
ramener à un problème à deux dimensions (sec- 
tions droites du cylindre). Dans de tels cas, aux- 
quels nous nous limiterons à partir d'ici, les équa- 
tions (5) et (4,a) sont à remplacer par : 


(5,a)‏ مر ربج ہے 
Ox? ay?‏ 
dV‏ 
4b‏ 29 
J (4,b)‏ 


désignant la dérivée du potentiel suivant la 


normale à la ligne équipotentielle, Quant aux con- 
ditions aux limites, elles s'ċnoncent comme plus 
haut, avec cette seule différence que le mot « sur- 
face » est à remplacer par le mot « ligne ». 


Lorsqu'il s'agit d'un problème à deux dimen- 
sions, la détermination de la répartition primaire 
se fait le plus souvent par la méthode des trans- 
formations conformes. 


Dans cette méthode, on considère deux variables 
complexes : 
Zee EE 


Z = EN 


Après avoir tracé dans le plan xy la configu- 
ration donnée, on essaie de la faire apparaître dans 


le plan X Y , grâce à une transformation appro- 
prices 


Z = f(z) (8) 


sous forme d'une configuration pour laquelle on 
connait la répartition primaire, par exemple sous 
la forme indiquée à la figure 5 (répartition uni- 
forme). 


Dès lors, partant des lignes équipotentielles et 
des lignes de courant de la répartition connue, on 
obtient les lignes homologues de la répartition cher- 
chée en utilisant la transformation inverse de (8) : 


zu E 9 (9) 


C'est ainsi que le problème d'une plaque plane 
(ici segment de droite), placée à l'intérieur d'une 
anode cylindrique de rayon infini, se ramène au 
cas de la répartition uniforme grâce à la transfor- 
mation. 

z = sinZ 


Le diagramme de champ que l'on obtient dans 
ce cas est représenté à la figure 6. La figure 7 
donne la courbe de répartition du courant sur une 
face de la plaque. Remarguons ici que cette courbe 
montre à quel point il est dépourvu de sens de 
rapporter les rendements cathodiques à la densité 
moyenne de courant ; celle-ci n'est égale à la den- 
sité locale qu'en deux points de la plaque (seg- 
ment). 


Bien que les méthodes analytiques, telles que la 
transformation conforme, permettent, en principe, 
de trouver, pour toute configuration donnée, la 
répartition primaire correspondante, leur appli- 
cation conduit à des calculs très compliqués dès 
que l'on s'ċcarte de formes simples. 


On comprend dès lors tout l'intérêt que présente 
l'emploi des méthodes analogiques : quelle que soit 


Du 


Champ autour d'une plaque plane avec l'anode à l'infini 
(R. RoussELoT) 


Fig. 6 
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Répartition du courant sur une face de la plaque plane 
(R. RoussEror) 


0 


la configuration à étudier, la répartition primaire 
est obtenue par un relevé direct. 


Parmi les procédés analogiques utilisés pour 
l'étude des champs plans, celui du papier conduc- 
teur est d'une mise en œuvre particulièrement 
commode. 


1.3.2 — Méthodes analogiques. Procédé du pa- 
pier conducteur (*) 


Dans ce procédé, l'ċlectrolvte est représenté par 
du papier conducteur. 


Pour obtenir le modèle analogique du système 
étudié, on découpe dans le papier, à une échelle 
donnée, la forme de ce système et l'on trace sur 
la découpe, à l'aide d'une peinture conductrice à 
base d'argent, des traits géométriquement sembla- 
bles aux électrodes du système. 


Lorsque celui-ci comporte (voir plus loin) des 
électrodes auxiliaires et des écrans isolants, les 
électrodes sont reproduites par des traits à la pein- 
ture conductrice et les écrans par des entailles pra- 
tiquées dans le papier conducteur (les parois de la 
cuve sont reproduites par les bords mêmes du 
papier découpé). 


Le modèle une fois terminé, on porte les électro- 
des aux potentiels voulus et l'on explore à l'aide 


(*) Cette méthode analogique a été systématiquement uti- 
lisée et développée par R. RoussEror (C.N.R.S. - Paris). 
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d'une sonde (simple pointe métallique à bout 
arrondi) le papier conducteur. 


Le montage électrique utilisé est celui schéma- 
tisé à la figure 8. 


Pour tracer une équipotentielle particulière, le 
curseur du potentiomètre E est placé dans la posi- 
tion correspondante. On déplace alors la sonde 
sur le papier jusqu'au moment où le galvanomètre G 
indique zéro ; la pointe de la sonde se trouve ainsi 
en un point de l'équipotentielle; on procède d'une 
maniére analogue pour obtenir d'autres points de 
l'équipotentielle. 

En pratique, l'objectif principal du galvano- 
chimiste est de connaître la distribution de la den- 
sité de courant à la surface de la cathode puisque 
l'épaisseur du dépôt en un point de celle-ci dépend 
de la densité de courant qui y existe. Comme la 
densité de courant est proportionnelle au champ 


électrique c'est-à-dire à [voir formule (4,b) ], 


le galvanochimiste est surtout intéressé par la 
connaissance du tracé d'une équipotentielle proche 
de la cathode. 


Pour passer à la densité de courant dans le bain 
d'électrolyse, il y a lieu de tenir compte du rap- 


Principe du tracé des équipotentielles sur un modèle en 
papier conducteur (modèle direct) (R. RousseLoT) 


Fig $ 


port de similitude géométrique et des rapports de 
tensions et de conductivités. 
Soient : 
I. : courant total dans le bain d'électrolyse (par 
unité d'épaisseur), 
L.: longueur caractéristique de la cuve d'élec- 
trolyse, 
ie: densité de courant dans le bain d'électrolyse, 
I, : courant total dans le modèle en papier, 


Lp: longueur caractéristique, homologue de Le, 
dans le modèle en papier, 

ip: densité de courant dans le modèle en papier 
(par unité de longueur), 

pp: résistivité de la couche conductrice, 

e: épaisseur de cette couche, 

AV, : différence de potentiel entre la cathode et 

une équipotentielle proche, 


An, : distance entre un point de la cathode et cette 
équipotentielle. 


Par un raisonnement élémentaire, on peut mon- 
trer que: 


Le 
i, = وا کا‎ )10( 
I, 
avec 
اھ‎ 
A 1 
ip TA (11) 
et 
L, یں‎ 
K — L (rapport de similitude 


géométrique) 
R, étant la résistance d’un carré de papier de 
coté a, on a: 
a Pp 
R, — Pp TOME 
ae e 


indépendant de a. (Pour le papier du type L, 
R, = 2.500 Q par carré). 


0 SC Inversion. 


Il est bien connu que le rċseau des lignes ċgui- 
potentielles d'une fonction harmonigue V(x,v) 
peut être considéré comme le réseau des lignes de 
force d'une autre fonction harmonique U(x , y); 
inversement, le réseau des lignes de force de 
V(x,y) s'identifie avec le réseau des équipoten- 


tielles de U(x , y). 


Pour réaliser une configuration de potentiel 
U(x , y), il suffit de construire un modèle dit 
inversé dans lequel, d'une part, les électrodes épou- 
sent le contour des parois isolantes du modèle 
initial et, d'autre part, les parois isolantes ċpousent 
le contour des électrodes. 


Utilisation d'une bi-sonde pour déterminer les densités de 
courant relatives sur un modèle inversé (R. RoussELOT) 


Fig. 9 


L'intérêt pratique du modèle inversé est qu'il 


, c'est-à-dire de 


V 
permet une mesure directe de d 

n 
la densité de courant au voisinage immédiat des 
électrodes. 


` 


Le moyen le plus simple consiste à utiliser une 
bi-sonde, c'est-à-dire un ensemble de deux sondes 
à écartement l faible et fixe (environ 5 mm) 
(fig. 9), reliée à un voltmètre à haute impédance. 
Si les pointes de la bi-sonde sont posées le long du 
contour d'une paroi isolante du modèle inversé 
(électrode du modèle direct) la différence de 
potentiel entre les pointes est une mesure de la 
densité de courant sur l'électrode du modèle direct. 


Si i, est le rapport : 
différence de potentiel entre les deux pointes de la bisonde 
différence de potentiel entre électrodes du modèle inversé 


on démontre que la densité de courant ie vaut : 


1.3.2.2. — Electrodes auxiliaires et écrans isolants. 


L'un des buts que recherche le galvanochimiste 
est de réaliser des dépôts d'épaisseur aussi uni- 
forme que possible sur des pièces de forme plus ou 
moins complexe, 
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La répartition de l'épaisseur est uniquement 
fonction de la répartition secondaire du courant 
car celle-ci s'établit dès les premiers instants de 
l'électrolyse et elle ne subit plus guère de modifi- 
cation ultérieure. Or, la polarisation, qui est res- 
ponsable de la répartition secondaire, est une gran- 
deur, en l'occurrence une différence de potentiel, 
qui se définit en chaque point de l'électrode et qui 
n'est fonction que de la densité de courant (secon- 
daire). 


On concoit donc qu'il soit possible d'amċliorer 
l'uniformité d'un dépôt si l'on peut agir sur les 
propriétés de polarisation de la solution. 


Effectivement, il existe aujourd'hui toute une 
gamme de substances qui remplissent cette fonc- 
tion ; il en sera question plus loin. 


Mais, d'autre part, il est bien clair que si la 
répartition primaire du courant était uniforme, il 
n'y aurait aucune raison pour que la répartition 
secondaire ne le soit pas; ainsi, le galvanochimiste 
a tout intérêt à réaliser une répartition de courant 
primaire, sur la cathode, aussi uniforme que pos- 


sible. 


On peut parfois améliorer cette répartition en 
plaçant au voisinage de la cathode, des écrans 
isolants ou des électrodes auxiliaires (reliées ou 
non à des sources de tension) : en effet, ces écrans 


et ces électrodes modifient les conditions aux limi- 
tes du potentiel V et, par là, la distribution de la 
dérivée normale de V le long de la cathode. 


À titre d'exemples, les figures 10, 11 et 12 illus- 
trent l'effet d'un écran isolant ou d'une électrode 
auxiliaire sur les réseaux d'équipotentielles et de 
lignes de courant dans le cas d'une cathode et 
d'une anode planes et parallèles. Ces figures ont 
été obtenues par R. ROUSSELOT à l'aide de modèles 
en papier conducteur. 


Il est à peine besoin de souligner que, dans la 
recherche d'une disposition favorable d'écrans et 
d'électrodes auxiliaires dans des structures cylin- 
driques, cette méthode analogique est d'une utilité 
absolument remarquable. 


I. 4. Détermination de la répartition secon- 
daire du courant. 


1.4.1 — Généralités. 
Le potentiel en un point P de l'électrolyte, 
proche d'une électrode, a pour expression : 
M(E) = va یت‎ E. — Es (7) 


dans laquelle V., potentiel de l'électrode, et E. , 
tension d'ċlectrode, sont des constantes ; par con- 


Cathode et anode planes et parallèles ; on observe la concentration des lignes de force au voisinage de l'extrémité de 
la cathode (R. RoussELor) 


Fig. 10 
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Cathode et anode identiques à celles de la fig. 10, avec un écran perpendiculaire à la cathode, placé près de l'extrémité 
de cette dernière ; on remarque une meilleure uniformité, d'après l'allure de la première équipotentielle près de la cathode 
(R. RousseLoT) 


Eig. li 


Cathode et anode identiques à celles des fig. 10 et 11, avec une électrode auxiliaire (anneau de garde) reliée à la 
cathode ; on note également une meilleure uniformité d'après l'allure de la premiere ċguipotentielle (R. RoussELor) 


7 12 
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tre, la polarisation e. est une quantité qui, ici, 
dépend de la densité de courant. 


Si l'on considère les courbes de la figure 3 et si 
l'on tient compte de la relation (4,b), on s'aperçoit 
que V(P) est une fonction non linéaire de —— 


dn 


au point P. 


Ainsi, la détermination analytique de la répar- 
tition secondaire du courant consiste, au point de 
vue mathématique, à rechercher une fonction har- 
monique satisfaisant à de telles conditions aux 
limites. C'est là un problème quasi insoluble sauf 
dans des cas extrêmement simples sans aucun 
intérêt pratique ou théorique. 


La méthode analogique au papier conducteur, 
elle-même, ne se préterait à une application directe 
que si Ton disposait de matériaux ou d'éléments de 
circuits électriques susceptibles de représenter les 
relations entre la polarisation et la densité de 
courant (*). 


I. 4.2. Méthode de R. ROUSSELOT 


R. ROUSSELOT a préconisé et développé un pro- 
cédé itératif basé néanmoins sur la méthode analo- 
gique au papier conducteur. 


En voici le principe. 


Il opère successivement sur plusieurs modèles 
en papier dont chacun correspond à une itération. 
Dans chaque modèle, sauf le premier, la différence 
de potentiel e., due à la polarisation cathodique, 
est représentée, à une certaine échelle, par une 
chute de tension e, due au passage du courant ip, 
proportionnel à i. , dans une résistance de valeur 
appropriée, qu'il s'agit chaque fois de déterminer. 
Cette résistance est réalisée de la manière suivante. 


Le long du contour correspondant à la cathode, 
R. ROUSSELOT prolonge le papier par une série de 
bandelettes de ce même papier, bandelettes qui 
sont rectilignes, de largeur uniforme A, isolées les 
unes des autres sur toute leur longueur. Les extré- 
mités libres de ces bandelettes sont alors réunies 
par un court-circuit sous forme d'une couche de 
peinture à l'argent. La figure 13 aide à comprendre 
cette disposition. 


Si e, est la chute de tension à réaliser dans une 
bandelette particulière tandis que i, est la densité 
de courant, la longueur l de cette bandelette est 
donnée par la formule : 


| = = (12) 


(*) A cet égard, l'emploi des varistances, ou résistances 
dépendant du courant, paraît être intéressant. 
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Modèle en papier conducteur avec analogie de la polarisation 
cathodique 


e 
Le rapport —— doit être constant et donne 
p 
l'échelle des tensions. 


Si V, est la somme de la chute de tension dans 
l'électrolyte et de la polarisation cathodique, et si 
V, est la chute de tension dans le modèle en 
papier, augmentée de la chute de tension dans la 
bandelette considérée, on doit donc avoir : 


CCE 7 
ee یت‎ (13) 


Dans la suite, cette relation sera utilisée dans 
les conditions suivantes : 


— V., et e. sont des données expérimentales, 


— V, est arbitraire ; 
on en déduira e, 


Il s'agit d'abord de mesurer V,. Dans ce but, 
on opère sur le modèle réel, c'est-à-dire dans la 
cuve électrolytique ; on place devant l'anode et très 
près de celle-ci une toile métallique de même métal 
que celui à déposer ; puis, on provoque l'électro- 
lyse jusqu'à l'établissement de la répartition secon- 
daire; on note le courant total I. (par unité 
d'épaisseur). 


Dans ces conditions, la différence de potentiel 
entre la toile métallique et la cathode n'est autre 
que V, car la polarisation anodique est éliminée 
de la mesure de même que les tensions d'électrode 
de la toile et du métal déposé sur la cathode 
(métaux identiques). 


Cela étant, on procède à la première itération. 


Celle-ci va conduire à une première valeur 
approchċe i,, de la densité de courant en chaque 
point de la cathode. 


Dans ce but, on utilise un modèle en papier 
identique à celui destiné au relevé de la distribution 
primaire de courant, ce qui revient à dire que la 
longueur des bandelettes est nulle dans ce premier 
modèle. 


La densité de courant i,, qui est mesurée est donc 
exactement celle correspondant à la distribution 
primaire ; par conséquent, ce relevé peut se faire 
sur un modèle inversé. Pour la détermination du 
courant i à partir du courant ipo , On utilise, non 
pas la formule (10) mais la formule équivalente : 


R, Ms 
K 114; 
A (14) 


dans laquelle K est encore le rapport de similitude 


leo = lpo 


p 


géométrique et pe la résistivité de l'électrolyte. 


La première itération est ainsi effectuée. 


Dans l'itération suivante, la valeur de ico , en tout 
point de la cathode, est utilisée pour déduire une 
première série de valeurs de la longueur des ban- 
delettes. 


A cette fin, la valeur de ico est portée en abscisse 
sur la courbe expérimentale donnant la polari- 
sation e. en fonction de i. (il s'agit donc d'une 
courbe analogue à celles de la figure 3) ; on obtient 
ainsi eo. 


La formule (13) donne alors une valeur appro- 
chée e,, de la chute de tension à réaliser dans une 
bandelette : 


€p1 = と co (15/2) 
Connaissant e et ip, on déduit de la for- 

mule (12) la valeur approchée 1, de la bandelette 

envisagée : 

لوت 


R, 158 


Dès lors, on peut passer à la réalisation du 
premier modèle en papier, muni des bandelettes de 
longueur lı. 


(12,a) 


= ہا 


La tension V, est alors appliquée à ce modèle, 
c'est-à-dire, plus précisément, entre l'ċlectrode de 
court-circuit des bandelettes et l'anode. 


On procède alors au relevé de la nouvelle den- 
sité de courant i,, le long du contour correspondant 
à la cathode, ce qui peut se faire par une mesure 
directe de la chute de tension dans chaque ban- 
delette; de i, on passe à i à l'aide de la 
relation (14). 


Ainsi se termine la deuxième itération. 


La suite du processus se calque sur celui qui 
vient d'être décrit, à savoir que de i,1, on passe à 
ei à l'aide de la courbe (e. , ie), puis à ep, par 
la formule : 


Epa = Er (13,b) 
CENT 


R. ROUSSELOT a remarqué qu'en pratique, le nom- 
bre d'itérations nécessaires dépasse rarement 
4 OU 


Bien que ce procédé soit quelque peu fastidieux, 
il possède au moins le mérite d'être le seul qui ait 
permis jusqu'ici de prévoir la distribution secon- 
daire du courant. 


I. 5. Rendement cathodique. 


S'il est bien vrai que la répartition du métal est 
uniquement fonction de la répartition secondaire du 
courant, il ne faudrait pas croire que la vitesse de 
formation du dépôt métallique en un point de la 
cathode soit simplement proportionnelle à la den- 
sité de courant en ce point. En fait, une fraction 
de cette densité est consacrée à la décharge d'hv- 
drogène ou à d'autres phénomènes secondaires ; 
cette fraction dépend de la valeur de la densité de 
courant secondaire. 

On définit donc le rendement cathodique ne 
comme étant le rapport de la densité de courant, 
qui, d'après la loi de FARADAY, contribuerait à la 
formation du dépôt, à la densité secondaire. 


Pour de faibles densités de courant, n. est tou- 
jours proche de 1 ou 100 % ; la loi de FARADAY 
est alors vérifée. Le rendement cathodique dimi- 
nue régulièrement lorsque la densité de courant 
croit. 


La figure 14 indique la dépendance de y. vis-à- 
vis de i, pour une solution de cyanure de cuivre. 
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Courbe de rendement cathodigue d'un bain de cvanure de 
cuivre en fonction de la densité de courant (R. RoussELoT) 


Fig. 14 
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Nous avons ainsi passé en revue les facteurs qui 
interviennent dans la répartition finale du métal 
et sommes maintenant en mesure de comprendre 
les moyens d'action physico-chimiques dont dis- 
pose le galvanochimiste pour améliorer l'uniformité 
des dépôts. 


I.6. Répartition du métal. 


Nous allons d’abord reprendre l'examen de 
l'influence de la polarisation et rechercher dans 
quel sens agit celle-ci au point de vue de l'unifor- 
mité du courant ; nous discuterons ensuite de l'effet 
du rendement cathodique. 


La chute de tension ohmique Vx — Vu, le long 
d'une ligne de courant dans le bain (fig. 15), est 
donnée par: 


NA E Vi 


| 


M _ 
۱ E dl 
N 

= 一 V。 一 (E, + €. — E.—ee) 
en prenant 


Va = 0 


05 


Les grandeurs Ve , E, E sont des constantes 
dans tout le système ; quant à e, , il peut générale- 
ment être considéré comme constant. 


La chute de tension Vx — Vu apparaît donc 
comme la somme d'un terme constant: — V, — 
(E, + e, 一 E.) et d'un terme variable: e.. 
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D'après l'allure générale des courbes de la 
figure 3, on voit que la valeur absolue de e. croît 
de façon monotone avec la densité de courant. 
Comme e. est négatif, la chute de tension Vy — Vu 
est d'autant plus faible que la densité de courant 
au point M est forte. 


D'autre part, d'après la loi d'OHM appliquée à un 
tube de courant aboutissant en M, la densité de 
courant est proportionnelle à Vx — Vu. Il en 
résulte que, si nous considérons dans la répartition 
primaire une région de la cathode où la densité 
de courant est forte, il y aura apparition, en cette 
région, d'une polarisation forte qui aura pour effet 
de diminuer le courant lorsqu'on laisse s'établir la 
répartition secondaire. Ainsi, on peut comprendre 
que, par suite de la polarisation, la répartition 
secondaire de courant soit plus uniforme que la 
répartition primaire et que cette action unifor- 
misante soit d'autant plus prononcée que e. varie 
plus rapidement avec i。, c'est-à-dire que la pente 
de la courbe (e. , i.) soit plus forte. 


Du raisonnement précédent, il résulte aussi que la 
polarisation est un agent d'uniformisation d'autant 
plus efficace que la valeur movenne de la tension 
Vx — Vu est faible car, dans ce cas, l'importance 
relative de e. est plus grande. Le fait d'augmenter 
la conductivité d'un électrolyte peut donc aussi 
avoir pour effet d'améliorer l'uniformité de la 
répartition secondaire. 


Comme il a été dit plus haut en 1.3.2.2, on dis- 
pose de substances: complexants, tensio-actifs,... 
qui exercent une influence sur la polarisation. Leur 
addition au bain d'électrolyte peut donc rendre de 
grands services. 


On peut aussi, localement, diminuer la polari- 
sation par une agitation mécanique du bain car 
celle-ci favorise la diffusion et le renouvellement 
de l'électrolyte. 


Enfin, il est aisé de concevoir que l’uniformite 
du dépót puisse encore étre améliorée en tirant 
profit de la diminution du rendement cathodique 
ne avec l'accroissement de ie; en effet, les zones 
à faible densité de courant bénéficient d'un meil- 
leur rendement cathodique. Par conséquent, le 
choix du courant le plus favorable ou de la den- 
sité de courant la plus propice, celle-ci étant 
influencée par l'interposition d'écrans ou d'élec- 
trodes auxiliaires, fournit au galvanochimiste un 
moyen de contrôle supplémentaire. 
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R.H. RoussELOT 


Publication : 


R.H. RoussELOT 
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CINQUANTIEME ANNIVERSAIRE DE LA M.B.L.E 
SEANCE ACADEMIQUE EN PRESENCE DE 
SA MAJESTE LE ROI LEOPOLD Ill 


VIJFTIGSTE VERJAARDAG VAN DE M.B.L.E ħ 
ACADEMISCHE ZITTING IN TEGENWOORDIGHEID. 
VAN ZIJNE MAJESTEIT KONING LEOPOLD Ill 


l'occasion du cinquantième anniversaire de sa 

hdation, la M.B.L.E a organisé, en la salle de la Madeleine 
Bruxelles, une séance académique solennelle qui a eu 
le 21 mars 1961. 

tte cérémonie fut honorée de la présence de 

| Majesté le Roi Léopold Ill. 

rmi l'assistance on remarquait un très grand nombre de 
rsonnalitċs appartenant aux communautés européennes, 
gouvernement et aux pouvoirs publics, aux mondes 
litique, militaire, scientifique, financier, industriel et 


mmercial. 
presse avait accordé son attention toute spéciale à 


‚enement. 


Une importante délégation du personnel était présente. 

M. Félix Leblanc, président du Conseil d'Administration, 
souhaita la bienvenue aux invités et les remercia de 
l'intérêt qu'ils portent à la MB.LE. 


M. Marcel Hublou, administrateur-délégué, 
remercié le Roi Léopold de son soutien actif à l'égard 
de ceux que préoccupe la recherche scientifique, fit un 
historique du développement de la M.B.L.E. 


après avoir 


Il insista sur l'ampleur des investissements de la MBLE 
dans le domaine de la recherche : 
firme a ses propres bureaux d'études, qui occupent 

280 personnes. Elle affecte à ce budget plus de 10 % de 
son chiffre d'affaires. 


depuis douze ans la 


M. Hublou souligna ensuite l'importance vitale que 
revêt à ses yeux l'humanisation du travail dans une 
entreprise de cette importance. 

Il attira également l'attention de l'auditoire sur le rôle 
important que jouent la sélection du personnel, sa 
ormation et son intégration au sein de l'entreprise. 
Succédant à M. Hublou, M. Charles Cammaert, délégué 
des Conseils d'Entreprise, se plut à reconnaître 
l'excellence du climat social qui règne parmi le personnel 


ontant ensuite à la tribune, M. Jacques Van Der Schueren, 

inistre des Affaires Economiques, dans un brillant 
exposé, cita les deux prérogatives remarquables que 
5ossede la M.B.L.E : un esprit de promotion de la recherche 
scientifique et technique, une politique sociale faite « du 
ens, du souci, du respect de l'homme 9. 

e Ministre décrivit les difficultés et les écueils qui 
juraient pu compromettre le succès d'une telle politique : 
lifficulté de trouver des chercheurs ayant les qualités 
ntellectuelles et morales requises, de les fondre en une 
squipe homogene sans écraser leur personnalité, d'intégrer 


tous les efforts pour exécuter un programme múrement | 
pensé ; problèmes de financement d'un matériel | 
scientifique de plus en plus coûteux et se démodant si | 
vite qu'on ne peut en garantir la rentabilité ; aléas des | 
fabrications expérimentales et de leurs transposition sur l 
plan industriel. C'est le mérite de la M.B.L.E d'avoir su | 
triompher de ces difficultés. | 
Parlant de la politique sociale, l'orateur dit notamment: | 
« La promotion sociale, ce sont d'abord de bons salaires 
»cest d'abord un maximum de sécurité d'existence — | 
» c'est-à-dire de prospérité. Le progrès social peut être un! 
» alibi, ou une excuse. Le premier devoir social d'un che 
» d'entreprise, c'est la gestion aussi rationnelle, aussi 

> heureuse que possible, de l'organisme qui lui est confié 
» Le second, c'est la rémunération aussi large que possiblé 
» des services qui lui sont prestés. » 

Le Ministre termina en décrivant la politique sociale de 
la MBLE comme reposant essentiellement sur le respec 
de la personnalité et sur l'exercice bien compris de 
l'autorité. 


Une réception très animée a suivi cette séance académique 


` gelegenheid van de vijftigste verjaardag van haar 
richting heeft de M.B.L.E, in de Madeleinezaal te Brussel, 
1 plechtige academische zitting ingericht, die plaats 
ep op 21 maart 1961. 


ze plechtigheid werd opgeluisterd door de aanwezigheid 
١ Zijne Majesteit Koning Leopold Ill. 


der de aanwezigen bevonden zich een groot aantal 
sonaliteiten, behorende tot de Europese 
neenschappen, de regering en de openbare diensten, 
politieke, militaire, wetenschappelijke, financiċle, 
ustriċle en commerciċle wereld. 

pers had zeer bijzondere aandacht aan de gebeurtenis 
vid. 

| sterke afvaardiging van het personeel was aanwezig. 
Heer Félix Leblanc, Voorzitter van de Raad van 
eer, heette de genodigden welkom en dankte hen 

r de belangstelling die zij betoonden voor de M.B.L.E. 
Heer Marcel Hublou, Afgevaardigde-Beheerder, 
etste, na Koning Leopold zijn dank betuigd te hebben 
r dezes actieve steun aan hen die zich bezighouden 


met wetenschappelijk onderzoek, de ontwikkelingsgang 
van de MBLE. 

Spreker legde de nadruk op de omvang van de investeringer 
van de MBLE op het gebied van het onderzoek : 
sedert twaalf jaar heeft de firma haar eigen studieburelen, 
waarin 280 personen werkzaam zijn. Aan dit budget 
besteedt zij meer dan 10 % van haar zakencijfer. 


De Heer Hublou beklemtoonde vervolgens het vitale 
belang van de humanisering van de arbeid in een 
onderneming van dergelijke omvang. 

Hij vestigde eveneens de aandacht van het gehoor op he 
grote belang van de selectie, de opleiding en de integratie 
van het personeel in de schoot van de onderneming. 

Na de Heer Hublou kwam de Heer Ch. Cammaert, 
afgevaardigde van de Bedrijfsraden aan het woord, die 
getuigenis kwam afleggen van het uitstekende sociaal 
klimaat dat heerst onder het personeel van de MB.LE. 
Daarna kwam de Heer Jacques Van Der Schueren, 
Minister van Economische Zaken op het spreekgestoelte 
en gaf een schitterende uiteenzetting over de twee grote 
voorrechten die de MBL E tot de hare heeft gemaakt : 


de geest van bevordering van het wetenschappelijk en 
technisch onderzoek en een sociale politiek die berust 
op «de zin, de bezorgdheid voor de eerbied voor de mens». 
De Minister schetste de moeilijkheden en de hindernissen 
die het succes van een dergelijke politiek hadden kunnen 
in gevaar brengen : moeilijkheid om vorsers te vinden, 
die de nodige intellectuele en morele eigenschappen 
bezitten, om ze te integreren in een ploeg zonder hun 
personaliteit te onderdrukken, om alle inspanningen samen 
te bundelen voor de uitvoering van een weloverwogen 
programma ; problemen van financiering van een 


wetenschappelijk materieel dat steeds duurder wordt en 
dat zodanig snel veroudert, dat de rentabiliteit er niet var 
kan gewaarborgd worden ; wisselvalligheden van de 


proeffabricages en van hun transpositie op industrieel plan. 
De M.B.L.E heeft de verdienste deze moeilijkheden te 
hebben kunnen overwinnen. 

Sprekende over de sociale politiek zegde spreker o.m. : 
« Sociale vooruitgang betekent in de eerste plaats goede 
» lonen en een maximum aan bestaanszekerheid, dat w 
» zeggen welvaart. De sociale vooruitgang kan een alibi 

» of een verontschuldiging zijn. De eerste sociale plicht 
» van een bedrijfsleider is een zo rationeel, een zo goed 
» mogelijk beheer van het organisme dat hem werd 

» toevertrouwd. De tweede is de zo goed mogelijke 

» bezoldiging van de diensten die hem worden geleverd. 
De Minister besloot zijn toespraak met een uiteenzetting 
over de sociale politiek van de MBLE, die hoofdzakeli 
berust op de eerbied voor de persoonlijkheid en op een 
goed begrepen uitoefening van het gezag. 

Na de academische zitting had een zeer levendige 
receptie plaats. 


décembre 1960): 
«On voit donc ainsi l 
» portée et la direction de 
» l'effort consenti par une 
» de nos grandes entrepris 
» ses électroniques qui a 
» rogné sur ses occupations 
» proprement commerciales le 
> temps et l'enthousiasme né- 
» cessaires à la création de ce 
» « Centre d'Information et de Ser- 
» vice pour l'Enseignement scienti- 
» fique et technique», qui entend 
> mettre, à prix coútant, ses réalisa- 
» tions et son matériel à la disposition 
» de notre enseignement national. 
» Comme elle le dit elle-même, cette 
» entreprise est une manière d'invitation 
» aux industriels, aux pédagogues et à tous 


l’industrie au tableau noir... » ceux qui veulent contribuer aux dévelop- 
» pement technique et scientifique, « seul 

La presse spécialisée a fait » capable d'orienter et d'harmoniser le déve- 
largement écho aux manifestations » loppement industriel en élargissant le cadre 
organisées récemment par » de leur vie ». Il y a, dans ce souci, une har- 

la M.B.L.E pour la présentation, » diesse qui doit être soulignée en même temps 
aux autorités de l'enseignement » qu'une affirmation de l'attention que portent les 
moyen et technique, du matériel » entreprises elles-mêmes à promouvoir un ensel- 
didactique qu'elle a développé » gnement qui constitue une véritable formation des 
pour faciliter l'enseignement » «hommes nouveaux» aux tâches que le progres 

de l'électronique. » apporte avec lui. Le but, en définitive, est d'empé- 

Des nombreux articles parus » cher qu'un fossé ne se creuse entre l'enseignement et 

à ce sujet, nous extrayons : » l'industrie, entre la formation et la pratique des métiers. 


» Cela coûte du temps et cela réclame des investissements 
» dont la rentabilité, certaine, est à échéance. 
» C'est-à-dire qu'il faut avoir foi en l'avenir pour entrepren- 
» dre cette œuvre et la mener à bien. C'est un autre honneur 
» de l'industrie belge que de s'en préoccuper. » 
@ de la revue « Toute la Radio » — édition belge — n° 253, sous 
la signature de M. H. D'Haese : 
« La M.B.L.E, dont le nom fait autorité en matière d'électronique, 
» s'est attachée à la réalisation d'un vaste programme de diffu- 
» sion d'informations techniques adaptées à tous les degrés de 
» connaissance. 
» L'effort réalisé par la M.B.L.E pour aider l'enseignement technique se 
» développe journellement ; la grande firme électronique a mis au point 
» un vaste programme de fabrication de matériel didactique, provisol- 
» rement limité à l'électronique et à la physique nucléaire. 

» Ce programme didactique a été mis au point avec la collaboration des 
» grandes industries électroniques et des représentants de tous les 

> secteurs de l'enseignement de plusieurs pays européens. 
» Il a été conçu pour faciliter la tâche de ceux auxquels on a 
» confié la formation des électroniciens et des physiciens 
» dont dépend l'avenir industriel et économique de la 

» Belgique. 
» Il réunit d'une part, la longue expérience et 
» les connaissances d'ingénieurs spécialisés 
» grâce auxquels la technique électro- 
» nique européenne s'est imposée 
» dans le monde entier, par son 
» standing élevé. » 


Bulletin van het V.B.N. 
(n° 31 van 10 decem- 
ber 1960): 


« 


» 


Zo ziet men dus het doel 

en de richting van deze 

insponning van een onzer 

grootste ondernemingen 

op electronisch gebied, die 

naast haar commerciéle oc- 

cupaties de tijd en het enthou- 

siasme heeft gevonden om een 

« Inlichtings- en dienstcentrum 

voor het wetenschappelijk en 
technisch onderwijs» op te rich- 

ten, dat tegen kostende prijs zijn 
verwezenlijkingen en zijn materieel 

ter beschikking van ons nationaal 
onderwijs wil stellen. Zoals zij zelf zegt, 
is dit initiatief een soort uitnodiging aan 
het adres van de industriélen, de op- 
voedkundigen en allen die willen bijdra- 


gen tot de technische en wetenschappelijke de industrie achter het schoolbord... 


ontwikkeling, «die alleen in staat is de 

indusiriéle vooruitgang te oriénteren en te 

harmoniseren door hun horizon te verruimen ». 

Er spreekt uit deze bezorgdheid niet alleen 

durf ; zij is tevens een bevestiging van de aan- 

dacht welke de ondernemingen zelf besteden aan 

het bevorderen von een onderwijs dat waarlijk 

« nieuwe mensen » kweekt voor de taken welke de 
vooruitgang medebrengt. Het uiteindelijke doel : ver- 
hinderen dat een kloof ontstaat tussen het onderwijs 

en de nijverheid, tussen de opleiding en de beroepsuit- 
oefening. Dat kost tijd en dat eist investeringen, die 
stellig — hoewel niet onmiddellijk — lonend zullen zijn. 
Dat wil dus zeggen dat men vertrouwen in de toekomst 
moet hebben om dit werk te ondernemen en tot een goed 
einde te brengen. Ook hier komt de Belgische nijverheid 
weer alle eer toe voor de wijze waarop zij zich ermede bezig 


houdt ». 


tijdschrift < Toute la Radio» — Belgische uitgave — n" 253, 
onder de handtekening van de h. H. D'Haese: 


« 


De M.B.L.E, waarvan de naam gezaghebbend is op het gebied 
van de elektronica, heeft zich tot taak gesteld een uitgebreid pro- 
gramma uit te werken van verspreiding van technische informatie, 
aangepast aon alle trappen van kennis. 


De gespecialiseerde pers heeft 
ruime aandacht gewijd aan de 
manifestaties, die onlangs door de 
M.B.L.E werden ingericht voor de 
voorstelling aan het middelbaar 
en technisch onderwijs, van het 
didactisch materieel dat zij heeft 
ontwikkeld voor de bevordering 
van het onderwijs in de elektronica. 
Uit de vele artikels die over dat 
onderwerp verschenen zijn, halen 
wij aan: 


» De inspanning van de M.B.L.E om het technisch onderwijs te helpen, 

» ontwikkelt zich dagelijks; de grote firma voor elektronica heeft een 
» uitgebreid fabricage programma uitgewerkt voor de produktie van 
» didactisch materieel, voorlopig beperkt tot de elektronica en de kernfysica. 
» Dit didactisch programma werd uitgewerkt in samenwerking met grote 

5 elektronische ondernemingen en vertegenwoordigers van al de sec- 

» toren van het onderwijs van verschillende Europese landen. 

» Het werd opgesteld om de taak te vergemakkelijken van hen, 

» aan wie de opleiding werd toevertrouwd van de elektro- 

» nici en fysici, van wie de industrieële en economische 

> toekomst van België afhangt. 

» Het verenigt enerzijds de lange ondervinding 

»en de kennis van gespecialiseerde inge- 

» nieurs, dank zij wie de Europese techniek 

»in zake elektronica, door haar hoge 

» standing, haar plaats veroverd 

» heeft in de gehele wereld. 9 


nouvelles extensions de la M.B.L.E 


La M.B.L.E crée une nouvelle 
unité de production à Roulers. 
Les études préliminaires concer- 
nant le marché du travail et l'im- 
plantation ont été faites par le 
«West-Vlaamse Economische 
Raad ». La fabrication de piéces 
pour l'électronique sera enta- 
mée cette année encore. Fin 1961, 
l'entreprise offrira déjà 300 em- 
plois. 


Le siège de Roulers sera la cin- 
quième unité de production de 
la M.B.L.E en Belgique; la qua- 
triċme, qui sera mise en marche 
sous peu, est située à Bruxelles, 
Boulevard Maurice Herbette. 


nieuwe uitbreiding van de mM.B.L.E 


De M.B.L.E roept een nieuwe 
produktie-eenheid in het leven te 
Roeselare. De voorbereidende 
studies betreffende de arbeids- 
markt en de implantatie werden 
gedaan door de « West-Vlaamse 
Economische Raad». Dit jaar 4 
nog wordt er de fabricage van 

onderdelen voor de elektronica 

aangevat. Einde 1961 zal de 

onderneming reeds werkgelegen- 

heid bieden aan 300 personen 


De zetel te Roeselare zal reeds 
de vijide produktie-eenheid van 
de M.B.L.E zijn in Belgie; de 
vierde die eerlang in bedrijf zal 
worden genomen, is gelegen te 
Brussel, Maurice Herbettelaan. 


ف 


25-30 MAART 1961 
ERNATIONAAL CENTRUM ROGIER 


IERPLAATS یں‎ BRUSSEL 


IXIV* Salon «INTERELECTRONIC» a eu lieu à Bru- 
es du 25 au 30 mars 1961. 


stand de la M.B.L.E, consacré à la pièce détachée, 
sentait quelques nouveautés qui se sont ajoutées 
on programme de vente. 


‚sont, entre autres: 


lun nouveau tube de sortie lignes, le PL 500, dont la 
particularité est une construction spéciale de l'anode 
qui réduit considérablement le courant secondaire 
de la plaque ; 


les tubes images rectangulaires de 23” et de 19”; 


tous les nouveaux accessoires pour tubes rec- 
tongulaires (transformateurs de sortie lignes et 
image, bobines de déflexion, contrôle de linéarité, 
Siet 


le nouveau sélecteur de canaux «MM» (Magic 
Memory), dans sa nouvelle exécution; 


un nouveau modèle de haut-parleur à encombre- 
ment réduit pourvu d'un blindage magnétique. 
Ce haut-parleur, spécial pour les équipements 
de TV, peut être placé 
n'importe où dans le récep- 
eur; 


les nouveaux combi-élé- 
ments, véritables briques 
de l'édifice électronique : 

des auto-transformateurs à 


rapport variable de toutes 
puissances ; 


quelques échantillons de 

l'impressionnante gamme 

des appareils de mesure Weston (U.S.A.) dont la 
M.B.L.E a la représentation exclusive pour la Bel- 
gique et le Grand-Duché de Luxembourg ; 


un petit échantillonnage du très vaste choix des 
tubes professionnels ADZAM. 


« Master Trainer » ou ensemble didactique à pan- 
rux amovibles y était également présenté et a 
enu tout le succès qu'il mérite. 


LE € RONIC (F.A.I.R.) 


Het XIVe Salon «INTERELECTRONIC » heeft te Brussel 


plaats gehad van 25 tot 30 maart. 


De stand van de M.B.L.E, gewijd aan de onderdelen, 
stelde enkele nieuwigheden voor, die haar verkoops- 


programma zijn komen aanvullen. 


Het betreft o.m.: 


De « Master Trainer », een didactisch toestel met uit- 
neembare panelen, werd er eveneens voorgesteld en 
oogstte er het succes, dat hij verdient. 


een nieuwe lijn-uitgangsbuis, de PL 500, waarvan 
de bijzondere eigenschap de speciale constructie is 
van de anode, die aanzienlijk de secundaire plaat- 
stroom beperkt ; 


rechthoekige beeldbuizen van 23” en 19"; 


al de nieuwe bijbehorende onderdelen voor rechi- 
hoekige beeldbuizen (beeld- en lijnuitgangstrans- 
formatoren, deflectiespoelen, lineariteitscontrole, 
CHA: 


de nieuwe kanalenkiezer «MM » (Magic Memory) 
in zijn nieuwe ultvoering ; 


een nieuw model luidspreker met beperkte afmetin- 
gen voorzien van een mognetische afscherming. 
Deze luidspreker, die speciaal voor TV-toestellen 
werd gebouwd, kan gelijk 
waar in de ontvanger wor- 
den opgesteld ; 


nieuwe combi-elementen, 
bouwstenen van het elek- 
tronica-gebouw ; 


autotranstormatoren met 
veranderlijke verhouding 
voor alle vermogens ; 


enkele stalen uit de indruk- 
wekkende keuze meettoe- 
stellen van Weston (U.S.A.) waarvan de M.BLE 
de exclusieve vertegenwoordiging heeft voor België 
en het Groothertogdom Luxemburg ; 


een beperkte keuze uit de grote verscheidenheid 
van professionele elktronenbuizen ADZAM. 


PHYSIQUE NUCLEAIRE 


Le département « Physique Nucléaire » de la M.B.L.E vient de mettre sur 
le marché une chambre d'ionisation portative à réponse logarithmique. 


Cet appareil est le plus récent de ceux mis au point par la M.B.L.E pour le 
contrôle du réacteur BR 2 à Mol. 


Le circuit, entièrement transistorisé, est alimenté par des piles standard 
de 1,5 volt ou par le réseau. La sensibilité de l'appareil permet de mesurer 
des intensités de radiution allant de 0,1 mR par heure à 100 mR par heure. 


KERNFYSICA 


Het departement « Kernivsica » van de M.B.LE heeft zopas een draagbare 
ionisatiekamer met logorithmische karakteristiek op de markt gebracht. 


Dit toestel is het nieuwste in de reeks die door de M.B.L.E werden ontwik- 
keld voor de controle van reactor BR 2 te Mol. 


De elektronische kring is volledig getransistoriseerd en wordt gevoed door 
standaard-elementen van 1,5 volt of door het net. De gevoeligheid van het 
toestel is zodanig dat stralingsintensiteiten gaande van 0,1 mR per uur tot 
100 mR per uur kunnen gemeten worden. 


CONGRES INTERNATIONAL DE SPECTROGRAPHIE A LYON 


Le Groupement pour l'Avancement des Méthodes Spectrographiques 
(G.A.M.S.) organise à Lyon, du 5 au 10 juin prochain, le c IX COLLOQUIUM 
SPECTROSCOPICUM INTERNATIONALE ». 

A cette occasion, le département Spectrographie de la M.B.L.E a préparé 
deux communications sur quelques applications récentes du contrôle spec- 
irographique, ainsi que sur un sujet de grande actualité, à savoir l'auto- 
matisation dans l'enregistrement et la transmission des résultats d'analyse 
spectrale. 

Au cours de la même période se tiendra une exposition spécialisée de 
matériel de spectrographie, où la M.B.L.E présentera sur un stand de 64 m' 
ses tout derniers équipements et plus particulièrement le Spectrovide et la 
console de mesure Metronic entièrement automatique avec impression des 
résultats d'analyse. 

Cette manifestation qui réunira à Lyon toutes les personnalités du monde 
des laboratoires aura, croyons-nous, une très grande répercussion tant par 
l'importance de son organisation que par la valeur et la qualité des travaux 
qui y seront présentés. 

La M.B.L.E, qui est restée à l'avant-plan de toutes les techniques d'avant- 
garde dans ce domaine, se devait dy participer très activement. 


INTERNATIONAAL CONGRES VOOR SPECTROGRAFIE TE LYON 


De « Groupement pour l'Avancement des Méthodes Spectrographiques 9 
(CAMS) richt van 5 tot 10 juni eerstkomend, te Lyon, het «IX COLLO- 
QUIUM SPECTROSCOPICUM INTERNATIONALE » in. 

Te dezer gelegenheid heeft het departement <Spectrogrofie > van de 
M.B.L.E twee mededelingen voorbereid over enige recente toepassingen 
van de spectrografische controle evenals over een onderwerp van grote 
actualiteit, te weten de automatie in de registrering en de overseining van 
de resultaten van de spectrale analyse. 

Tijdens dezelfde periode zal een gespecialiseerde tentoonstelling worden 
gehouden van spectrografisch materieel, waar de M.B.L.E in een stand met 
een oppervlakte van 64 m' zijn laatste uitrustingen zal voorstellen en in 
het bijzonder de « Spectrovide » en de volledig automatische meetconsole 
« Metronic » die de analvseresultaten afdrukt. 

Deze monifestatie, die te Lvon alle personaliteiten uit de wereld van het 
laboratorium zal bijeenbrengen, zal, denken wij, een grote weerklank vinden, 
zowel door de omvang van haar organisatie als door de waarde en de 
kwaliteit van de werken die er zullen voorgedragen worden. 

De M.B.L.E, die op de voorposten is gebleven van alle avant-gardetechnieken 
op dat gebied, achtte het zich tot plicht daaraan zeer actief deel te nemen. 


De nederlandstalige lezers zullen hier- 
onder een bondige samenvatting vinden 
van de artikelen van dit nummer. 
Indien ze nadere inlichtingen wensen 
over de behandelde onderwerpen, wor- 
den ze verzocht zich in verbinding te 
stellen met « Manufacture Belge de 
Lampes et de Matériel Electronique 
N.V.» Twee Stationsstraat, 80, Brussel, 


wetenschappelijke en technische publicatie 


5de jaargang 4de deel 


Het praktisch gebruik van de Matrixrekening 
[Iste deel) 


door A. BRUAUX 


Na gewezen te hebben op de fundamentele 
rol, die de matrixrekening speelt in de 
moderne Fysica toont de schrijver het 
belang aan van het gebruik van de metho- 
den der matrix-algebra voor het oplossen 
van problemen, die de ingenieur en de 
fysicus in de praktijk te behandelen krijgen. 
Daartoe brengt hij de grondbeginselen van 
deze algebra in herinnering en bepaalt dan 
het karakter van operator, dat men kan 
toekennen aan de correspondentie-matrices. 
Ten slotte werkt hij twee zeer verschillende 
voorbeelden uit en toont aldus aan hoe 
deze matrix-methodes kunnen gebruikt 
worden en hoe soepel ze in het gebruik 


zijn. 
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Beginselen van de Meting der Stralingen 


(3de deel) 


door A. LAGASSE 


In deze nota zet-de auteur in grote trekken 
de weerslag uiteen van het wisselvallig 
karakter van de verschijnselen die optreden 
bij de meting van stralingen, op de uitslag 


van deze metingen. 


Hij onderzoekt eerst de oorsprong van dit 


wisselvallig karakter. 


Vervolgens formuleert hij de wet van de 


spreiding der probabiliteiten met betrekking 
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tot het aantal geregistreerde pulsen. In dit 
opzicht herinnert hij aan de grondbegin- 
selen van de wiskundige statistiek, die de 
gebruikte praktische methodes rechtvaardi- 
gen en toelaten de statistische fout te 
schatten bij een telling. Een korte ontleding 
van de invloed van de eigen beweging van 
de detectie-installatie op de nauwkeurigheid 
der metingen en enkele beschouwingen over 
de correcties te wijten aan de resolutie- 
tijd van deze installatie besluiten deze 
studie. 


ODE ESA 


Enkele uitzichten van de techniek der Gal- 


vanostegie (Iste deel) 


door J. DE RUDDER 


In deze eerste nota onderzoekt de auteur 
de verscheidene fundamentele grootheden 
zoals het diffusiepotentieel, de elektrode- 
spanning, de polarisatie. Hij brengt ver- 
volgens de wetten in herinnering die de 
verdeling van de stroom in de elektrolyten 


beheersen. 


Na een verwijzing naar de analytische me- 
thodes behandelt hij in het bijzonder de 
analogische methodes, die toelaten de ver- 


deling van deze stroom te voorzien. 


Ten slotte beschrijft de auteur de werking 
van de verschillende fysische en scheikun- 
dige faktoren die van aard zijn om de 
metaalneerslag te kunnen verbeteren uit 


het oogpunt van de eenvormigheid. 


WED GERO IES DAE 


REVUE NK KE KA KS 


English readers will find below a brief 
summary of the principal articles in 
this review. If they wish to have further 
information regarding the questions 
dealt with, they are requested to be 
good enough to communicate direct 
with the Manufacture Belge de Lam- 
pes et de Matériel Electronique, S.A., 
80, rue des Deux-Gares, Brussels. 


Practical applications of matrix calculus 
(Part 1) 


by A. BRUAUX 


Emphasizing the fundamental function being 
performed by matrix calculus in modern 
physics, the author stresses the interest 
which is derived from the use of matrix 
algebrical methods for the resolution of 
problems which are practically encountered 
by the engineer and physicist. 


In view of this, the author first reminds us 
the elements of such algebra and gives 
evidence of the characteristics of operator 
which can be conferred to correspondence 
matrix. 


Finally, he operates on two very different 
examples which illustrate the application of 
matricial methods and show their flexibility. 
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Principles of radiation measurements 
(Part 3) 


by A. LAGASSE 


In the present note, the author shows, in 
its broad lines, the incidence of the alea- 
tory characteristics of the phenomenons 
influencing the results of radiation measu- 
rements. 


He first analyses the origins of such alea- 
tory characteristics. 


He thereafter formulates the law which 
governs the probability distribution in relation 
to the number of recorded impulsions. He 
reminds the reader of the fundamental 
elements of mathematical statistics which 
justify the practical methods thus used and 
also the evaluation of the statistical error 
made on a numeration. 


A brief analysis of the influence of self- 
motion of the detection equipment on the 
precision of measurement and a few consi- 
derations on the correction due to the reso- 
lution time of such equipment terminate the 
study. 


Some aspects of electro-deposition techni- 


ques (Part |) 


by J. DE RUDDER 


The author reviews the fundamental units 
such as of diffusion potential, electrode 
tension and polarization. 


He reminds the reader of the laws which 
regulate the repartition of current in elec- 
trolytes. $ 


Continuing with a brief reference to analy- 
tical methods, the author more thoroughly 
examines the analogical methods which 
allow a prevision of such current distribution. 


The author also shows the modes of action 
of various physical and chemical factors 
which affect the improvement om metal 
deposition from the point of view of uni- 
formity. 


1961 


News in brief 


its fiftieth anniversary by a solemn acade- 
mic cession. 


of His Majesty King Leopold Ill. e 
Addresses were delivered by Mr. F.-LE- 
BLANC, Chairman, Mr. M. ۱٣ 
Managing Director and Mr. CAMMAERT, 
representing the personnel. 
Mr. J. VAN DER SCHUEREN, the belgian 
Minister of National Economy, made ê 
brillant speech on the main traits of the 
M.B.L.E : a spirit of promotion of scientific 
research and technical development cou- 
pled with a policy of social promotion 
based on a constructive conception of 
human relations. A 
The Press issued extensive reports on the 
event. : 


2. M.B.L.E has set up a new factory at 
Roulers, where the manufacture of elec- 
tronic components will be started this year. 
The Roulers factory will be the fifth pros 
duction unit of M.B.L.E (the fourth one 
will soon be operating in Brussels). 


3. The fourteenth « Interelectronic > exhi- 
bition took place in Brussels from March 
25th till March 30 th 1961. 

At the M.B.L.E stand, mostly devoted to 
electronic components, the new items just 
added to the sales program were shown to 
the specialist visitors. ; 


4. The Nuclear Physics Department of 
M.B.L.E is now marketing a new portable 
ionisation chamber with logarithmic response. | 


5. The «Groupement pour l'Avancement. 
des Méthodes Spectrogranhiques» (G.A. 
M.S.) is organising at Lyon, from June 5th 
till June lOth next, the «IX Colloquium — 
Spectroscopicum Internationale ». y 
The Spectrographic Department of M.B.L.E T 
will present two communications on spec- | 
trographic control and automation in the © 
recording and transmission of spectral ana- _ 
lysis. | 
M.B.L.E will take part, during the same — 
period, in an exhibition of spectrographic 
equipment where, among other equipments, 
the ««Spectrovide » and the Metronic mea- 
surement console will be displayed. 


